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R87-957333-1 


The  Unsteady  Flow  in  the  Far  Field  of  an  Isolated  Blade  Row 


SUMMARY 

This  report  describes  a  theoretical  investigation  of  the  unsteady  flow  in  the  far  field 
of  an  isolated  two-dimensional  blade  row.  In  particular,  analytic  far-field  solutions 
are  derived  for  cascades  operating  at  subsonic,  transonic  and  supersonic  inlet  and  exit 
Mach  numbers.  These  solutions  describe  the  velocity  potential  fluctuations  associated 
with  irrotational  pressure  (acoustic)  disturbances,  rotational  velocity  disturbances  and 
the  vorticity  shed  from  blade  trailing  edges  and  convected  along  blade  wakes.  They 
thus  provide  the  inlet  and  exit  boundary-condition  information  needed  for  a  proper 
numerical  resolution  of  the  unsteady  flow  through  the  cascade.  As  part  of  this  investi¬ 
gation  a  convenient  method  for  classifying  unsteady  excitations  based  on  the  acoustic 
response  in  the  far  field  is  developed  and  analytical  results  are  presented  to  illustrate 
subsonic  and  supersonic  far-field  acoustic  response  behavior.  Numerical  calculations 
have  also  been  carried  out  for  unsteady  flows  through  subsonic  compressor-type  cas¬ 
cades  to  demonstrate  the  results  of  matching  analytic  far-field  to  numerical  near-field 
solutions  and,  to  some  extent,  the  impact  of  unsteady  far-field  behavior  on  the  aero¬ 
dynamic  response  at  a  vibrating  blade  surface. 
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1.  INTRODUCTION 


There  is  an  ongoing  need  to  develop  more  complete  unsteady  aerodynamic  pre¬ 
diction  methods  for  aeroelastic  and  aeroacoustic  applications.  Currently,  linearized 
unsteady  aerodynamic  analyses  in  which  unsteady  fluctuations  are  regarded  as  small- 
amplitude  perturbations  of  a  nonuniform  isentropic  and  irrotational  mean  or  steady 
flow  are  receiving  considerable  attention  (for  recent  reviews  see  Atassi  1987  and  Ver- 
don  1987a, b),  particularly  for  turbomachinery  aeroelastic  applications  (Atassi  &;  Akai 
1980;  Caruthers  1981;  Whitehead  1982;  Verdon  &:  Caspar  1982,  1984).  In  this  situation 
the  unsteady  flow  is  excited  by  prescribed  small-amplitude  temporally  and  spatially  pe¬ 
riodic  structural  and  external  aerodynamic  disturbances,  and  the  complex  amplitudes 
of  the  dependent  linearized  unsteady  flow  variables  are  governed  by  a  system  of  lin¬ 
ear  time-independent  equations  which  contain  variable  coefficients  that  depend  on  the 
underlying  mean  flow. 

As  a  consequence  of  this  linearization,  boundary  conditions  at  moving  solid  surfaces 
and  jump  conditions  at  moving  wakes  and  shocks  can  be  imposed  at  the  mean  positions 
of  such  surfaces.  Also,  for  turbomachinery  applications  the  required  solution  domain 
can  be  limited  to  a  single  extended  blade-passage  region.  Moreover,  analytical  far-field 
solutions  can  be  determined  and  matched  to  numerical  near-field  solutions;  hence,  this 
solution  domain  can  be  restricted  further  to  a  single  blade  passage  region  of  finite  extent 
in  the  axial  direction.  These  features  lead  to  economical  and  therefore  technologically 
useful  unsteady  aerodynamic  response  predictions. 

In  this  report  analytic  solutions  for  the  unsteady  velocity  potential  response  far  up¬ 
stream  and  far  downstream  of  an  isolated  two-dimensional  blade  row  are  determined. 
We  consider  the  general  unsteady  cascade  problem  in  which  the  unsteady  flow  is  pro¬ 
duced  by  structural  (i.e.,  blade  motions)  and/or  external  aerodynamic  excitations  (i.e., 
incident  acoustic,  entropic  or  vortical  gusts)  and  provide  solutions  for  compressor-  and 
turbine-type  cascades  operating  under  subsonic,  transonic  and  supersonic  inlet  and 
exit  conditions.  Fourier  series  solutions  will  be  determined  for  the  continuous  velocity 
potential  fluctuations  associated  with  the  acoustic  response  in  the  far  upstream  and  far 
downstream  regions  as  well  as  for  the  continuous  potential  fluctuations  in  the  far  down¬ 
stream  region  associated  with  the  rotational  velocity  disturbances  that  are  convected 
through  the  blade  row.  A  closed  form  solution  will  be  determined  for  the  discontinuous 
potential  fluctuation  associated  with  the  concentrated  vorticity  convected  along  blade 
wakes.  The  specific  form  of  the  solution  for  the  potential  due  to  acoustic  response 
fluctuations  depends  upon  whether  the  free-stream  Mach  number,  M,  is  less  than  or 
greater  than  one,  the  free-stream  flow  angle,  fi,  is  positive  or  negative,  the  axial  Mach 
number,  M  cos  H,  is  less  than  or  greater  than  one  and  the  temporal  frequency,  a of  the 
unsteady  motion  is  positive  or  negative.  A  single  solution  holds  for  the  potential  fluctu- 
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ation  associated  with  rotational  velocity  fluctuations  and  concentrated  wake  vorticity 
regardless  of  the  sign  of  M  —  1,  fl,  M  cos  ft  —  1  or  w. 

The  complete  solution  for  the  velocity  potential  in  the  far  field  provides  the  inlet 
and  exit  information  needed  for  the  linearized  unsteady  aerodynamic  analyses  intended 
for  blade  aeroelastic  and  aeroacoustic  design  predictions.  Moreover,  a  knowledge  of  the 
unsteady  aerodynamic  response  in  the  far  field  is  a  prerequisite  to  guide  the  choice  of 
computational  grid  spacings  needed  to  resolve  accurately  the  unsteady  flow  through 
the  blade  row,  particularly  for  cascades  operating  at  high  subsonic  or  supersonic  inlet 
and  exit  Mach  numbers  or  subjected  to  high-frequency  unsteady  excitations.  Such 
information  should  also  be  useful  for  constructing  time-accurate  numerical  simulations 
(e.g.,  see  Rai  1985;  Fransson  &  Pandolfi  1986;  Giles  1987)  of  nonlinear  unsteady  flows 
through  turbomachine  cascades. 

Once  the  complete  far-field  solution  is  given  below,  the  far-field  acoustic  response  is 
studied  leading  to  a  useful  classification  of  unsteady  excitations.  In  addition,  the  results 
of  parametric  studies  are  presented  to  illustrate  the  acoustic  response  behavior  in  the 
far  field  at  various  free-stream  Mach  numbers  and  flow  angles.  Finally,  we  present 
numerical  solutions  for  unsteady  flows  through  subsonic  compressor-type  cascades  to 
demonstrate  the  results  that  can  be  achieved  by  matching  analytic  far-field  solutions 
to  numerical  near-field  solutions  and  to  illustrate  the  effects  of  the  far-field  acoustic 
response  on  the  aerodynamic  response  (i.e.,  surface  pressures  and  unsteady  airloads) 
at  a  vibrating  blade  surface. 
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2.  PHYSICAL  PROBLEM  AND  GOVERNING  EQUATIONS 

We  consider  the  time-dependent  two-dimensional  adiabatic  flow,  with  negligible 
body  forces,  of  an  inviscid  non-heat-conducting  perfect  gas  through  a  cascade  of  airfoils 
or  blades  such  as  the  one  shown  in  figure  1.  All  physical  quantities  are  dimensionless. 
Lengths  have  been  scaled  with  respect  to  blade  chord,  time  with  respect  to  the  ratio 
of  blade  chord  to  the  upstream  free-stream  speed,  density  and  velocity  with  respect  to 
the  upstream  free-stream  density  and  velocity,  respectively,  and  pressure  with  respect 
to  the  product  of  the  upstream  free-stream  density  and  the  square  of  the  upstream 
free-stream  speed.  The  mean  or  steady  state  positions  of  the  blade  chord  lines  coincide 
with  the  line  segments  77  =  £  tan©  +  mG,  0  <  £  <  cos©,  where  £  and  77  are  Carte¬ 
sian  axes  attached  to  the  blade  row  and  pointing  in  the  axial-flow  and  the  cascade 
“circumferential”  directions,  respectively,  m  =  0, 1,2, ...  is  a  blade  number  index,  0  is 
the  cascade  stagger  angle  and  G  is  the  cascade  gap  vector  which  is  directed  along  the 
77-axis  with  magnitude  equal  to  the  blade  spacing. 

In  the  absence  of  unsteady  excitation  the  flow  beyond  some  finite  distance  upstream 
(say  £  <  £_)  and  downstream  (£  >  £+)  from  the  blade  row  is  assumed  to  be  at  the  most 
a  small  irrotational  steady  perturbation  from  a  uniform  free  stream.  In  addition,  the 
unsteady  flow  remains  attached  to  the  blade  surfaces  and  therefore,  thin  vortex  sheets 
or  unsteady  wakes  emanate  from  the  blade  trailing  edges  and  extend  downstream. 
Finally,  any  shocks  that  might  occur  are  assumed  to  be  of  weak  to  moderate  strength 
and  have  small  curvature.  Thus,  changes  in  the  entropy  and  vorticity  of  a  fluid  particle 
as  it  passes  through  a  shock  are  regarded  as  negligible. 

The  time-dependent  or  unsteady  fluctuations  in  the  flow  arise  from  one  or  more 
of  the  following  sources:  blade  motions,  upstream  and/or  downstream  acoustic  distur¬ 
bances  which  carry  energy  toward  the  blade  row,  and  upstream  entropic  and  vortical 
disturbances  which  are  convected  through  the  blade  row.  These  excitations  are  assumed 
to  be  of  small  amplitude  and  periodic  in  time.  The  external  aerodynamic  excitations  are 
also  spatially  periodic,  while  the  blade  motion  is  periodic  in  the  cascade-  or  77-direction. 
For  example,  we  consider  blade  motions  of  the  form 

R  ( X  +777  G,  £)  =  Re{  r  (X)  exp[t(u;£  +  ma)],  X  on  B  .  (2.1) 

In  (2.1)  R  measures  the  displacement  of  a  point  on  a  moving  blade  surface  relative 
to  its  mean  or  steady-state  position,  X  is  a  position  vector,  t  is  time,  r  is  a  complex 
displacement-amplitude  vector,  u)  is  the  frequency  of  the  blade  motion,  a  is  the  phase 
angle  between  the  motions  of  adjacent  blades,  Re{  }  denotes  the  real  part  of  {  }  and 
B  denotes  the  reference  (777  =  0)  blade  surface.  Incident  disturbances  are  of  the  form 

Pi(X,t)  =  i2e{p/iT00 exp[t(/cT00  •  X  +  wf)],  f  ,  (2.2) 

where  pi^0 o  and  /cTOO  are  the  amplitude  and  wave  number,  respectively,  of  an  incident 
pressure  fluctuation,  pj(X,t),  coming  from  far  upstream  (  —  00)  or  far  downstream  (00) 
and  carrying  energy  towards  the  blade  row. 
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Note  that  the  interblade  phase  angle,  a,  of  an  incident  disturbance  is  ktoo  G.  Also, 
the  temporal  frequency  and  wave  number  of  an  incident  entropic  or  vortical  disturbance 
are  related  by  c a  =  — •  V-oo  where  F_oo  is  the  uniform  relative  inlet  velocity,  but 
a  more  complicated  relationship  exists  between  u  and  for  an  incident  pressure 

disturbance  as  will  become  clear  from  the  discussion  below. 

For  small-amplitude  unsteady  excitations,  i.e.,  |r|,  |p/(Too|,  etc.  ~  0(e)  <C  1,  the 
time-dependent  flow  can  be  regarded  as  a  small  perturbation  about  an  underlying 
mean  flow.  In  addition,  because  of  our  assumptions  regarding  shocks  and  the  flow  far 
upstream  of  the  blade  row,  the  underlying  mean  or  steady  flow,  which  is  assumed  to 
be  known  in  the  present  study,  will  be  isentropic  and  irrotational.  Thus  V  =  and 

=  -(I"'-''1  =  b 


=  1-  h  2 


2  +  h  -  1  )Mlm 

2  +  (7  -  1)M2 


(2.3) 


where  F,$,  p,  M,  P  and  A  are  the  local  steady-flow  velocity,  velocity  potential,  density, 
Mach  number,  pressure  and  speed  of  sound  propagation,  respectively,  and  7  is  the 
specific  heat  ratio  of  the  fluid.  Finally,  because  of  the  assumed  form  of  the  unsteady 
excitations,  the  first-order  flow  properties  must  be  harmonic  in  time,  e.g.,  p{X,t)  = 
Re{p(X)e'wt} ,  and  they  must  satisfy  the  blade-to-blade  periodicity  condition,  i.e., 


p(X  +m  G)  =  p(X)eima  . 


(2.4) 


Thus  the  complex  amplitudes  of  the  linearized  unsteady  flow  variables  can  be  deter¬ 
mined  by  solving  a  time-independent  set  of  governing  equations,  and  numerical  solu¬ 
tions  to  these  equations  are  required  only  over  a  single  extended  blade-passage  region. 
The  field  equations  governing  the  first-order  unsteady  flow  can  be  written  as  (Gold¬ 


stein  1978,  1979) 


(2.5) 


A-iv"  -sV»/2)  +  IP*  -sV$/2)  •  V|V$  =  0  ,  (2.6) 

and  _  _ 

A(A-j^)-r‘v-(w)  =  r‘v.(?rs) .  (2.7) 


Here  s,  vR  and  (f)  are  the  complex  amplitudes  of  the  first-order  entropy,  rotational 
velocity  and  velocity  potential,  respectively,  and  D / Dt  =  iu>  +  V  •  V  Is  a  mean  flow 
convective  derivative  operator.  The  complex  amplitude  of  the  first-order  velocity  is 
given  by  v  =  +  vR.  The  rotational  velocity  is  divergence  free  far  upstream  of  the 

blade  row,  and  the  unsteady  pressure  fluctuation  is  related  directly  to  the  velocity 
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potential  by  p  =  —pV<fi/Dt.  In  general  we  require  a  solution  to  the  foregoing  system  of 
equations  which  is  subject  to  the  condition  of  flow  tangency  at  blade  surfaces  (attached 
flow),  appropriate  jump  conditions  across  shocks  and  blade  wakes,  and  the  requirement 
that  acoustic  response  disturbances  either  attenuate  or  propagate  away  from  or  parallel 
to  the  blade  row  in  the  far  field  (Yerdon  1987a, b).  As  a  consequence  of  the  small 
unsteady  disturbance  approximation,  surface  conditions  can  be  imposed  at  the  mean 
positions  of  the  blade,  shock  and  “wake”,  i.e.,  the  downstream  stagnation  streamlines, 
surfaces  with  mean  shock  and  wake  locations  being  determined  from  the  nonlinear 
steady  solution. 

The  system  of  linearized  unsteady  equations  can  be  solved  sequentially  to  deter¬ 
mine  the  fluctuations  in  entropy,  rotational  velocity  and  velocity  potential  throughout 
an  extended  blade  passage  solution  domain.  The  entropic  and  rotational  velocity  fluc¬ 
tuations  are  convected  by  the  mean  flow  and  therefore  can  be  determined  in  terms  of 
prescribed  upstream  (i.e.,  on  £  =  £_)  entropic  and  rotational  velocity  distributions. 
Information  on  the  velocity  potential  fluctuation  is  required  in  both  the  far  upstream 
and  far  downstream  regions  of  the  flow.  Moreover,  this  information  cannot  be  provided 
a  priori  as  it  depends  on  the  response  of  the  blading  to  an  imposed  unsteady  excitation. 
However,  analytic  solutions  for  the  potential  in  the  far  field  can  be  determined.  These 
can  be  matched  to  a  (numerical)  near-field  solution  and  thereby  serve  to  complete  the 
specification  of  the  linearized  unsteady  boundary-value  problem.  This  feature  also  per¬ 
mits  the  numerical  solution  domain  to  be  limited  to  a  single  extended  blade  passage 
region  of  finite  extent  (i.e.,  £_  <  £  <  £+)  in  the  axial-flow  direction. 

2.1  Governing  Equations  in  the  Far  Field 

In  the  far  upstream  and  far  downstream  regions  the  steady  background  flow  is  at 
most  a  small  perturbation  (of  0(e))  from  a  uniform  free  stream.  Therefore  it  follows 
from  (2.7)  that  far  from  the  blade  row  and  to  within  first-order  in  e,  the  unsteady  po¬ 
tential  is  governed  by  the  constant-coefficient  differential  equation  of  classical  linearized 
unsteady  aerodynamic  theory  (see  Miles  1959;  Ashley  &  Landahl  1961);  i.e., 

=  (i  fe.  (2.8) 

It  also  follows  from  (2.6)  that  the  rotational  velocity  is  governed  by 

~ET  (*w+  '^)  vR=  0  •  (2-9) 

Note  that  under  our  present  assumptions  \/  •  vR  =  0  m  the  far  upstream  region  and 
therefore  the  far  upstream  velocity  potential  fluctuations  are  caused  by  acoustic  dis¬ 
turbances  only. 

We  will  require  analytic  solutions  to  (2.8)  for  periodic  (in  t])  but  otherwise  arbi¬ 
trary  potential  distributions  at  the  axial  locations  £  =  ^  and,  as  a  consequence,  to 


5 


(2.9)  for  periodic  but  otherwise  arbitrary  rotational  velocity  distributions  on  £  =  £+. 
These  solutions  must  satisfy  the  blade-to-blade  periodicity  condition  (2.4)  and  also  the 
conditions  of  continuity  of  pressure  and  normal  velocity  across  blade  wakes;  i.e., 

=  0  or  |<£]  =  exp\im(e  -  Tm/V «,)],  X  on  Wm,  £  >  £+  ,  (2.10) 

and 

[V0]-  n=  n,  X  on  Wm,  £  >  £+  ,  (2.11) 

respectively.  In  addition,  the  potential  <f>  must  obey  the  condition  that  acoustic  response 
waves  either  attenuate  or  propagate  away  from  or  parallel  to  the  blade  row  in  the  far 
field.  In  (2.10)  and  (2.11)  the  symbol  [  ]  refers  to  the  difference  (upper  —  lower)  in 
a  flow  quantity  across  a  wake,  |  J pej  is  the  value  of  [  J  at  the  point  of  intersection 
(£+,7y+)  of  the  reference  (m  =  0)  wake  and  the  axial  line  £  =  £+  and  Tm  is  a  coordinate 
measuring  distance  in  the  direction  of  V, downstream  from  the  line  £  =  £+.  Note  that 
to  within  first-order  in  e,  the  mean  wake  positions,  Wm,  far  downstream  of  the  blade  row 
can  be  regarded  as  an  array  of  straight  lines  parallel  to  the  free-stream  flow  direction. 
We  require  such  solutions  for  subsonic  ( MTOO  <  1)  and  supersonic  (A/Tco  >  1)  relative 
inlet  and  exit  conditions.  However  in  the  supersonic  case,  we  are  interested  primarily 
in  flows  having  subsonic  axial  velocity  component,  i.e.,  MToocosfITOO  <  1  since  this  is 
the  situation  of  practical  interest  for  axial-flow  turbomachines. 
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3.  UNSTEADY  FAR-FIELD  SOLUTIONS 


3.1  Preliminaries 

To  determine  solutions  for  the  unsteady  potential  in  the  far  field  it  is  convenient  to 
set 

*(*)  =  4,p(X)  +  ^(X)l/({  -  (+)  forfS^,  (3.1) 

where  cj>p(X)  is  a  continuous  function  which  accounts  for  the  fluctuations  in  the  ve¬ 
locity  potential  due  to  irrotational  pressure  or  acoustic  disturbances,  <pR{X)  accounts 
for  the  fluctuations  in  the  potential  in  the  far  downstream  region  due  to  rotational 
velocity  disturbances  (including  those  due  to  the  vorticity  shed  at  blade  trailing  edges 
and  convected  along  the  blade  wakes)  and  U(x)  is  the  unit  step  function.  The  latter 
disturbances  are  convected  as  vorticity  waves  by  the  underlying  uniform  flow  and  have 
no  associated  pressure  or  density  fluctuations  (Smith  1971).  Therefore, 

~Dt4’R=0  (3'2) 

and  it  follows  from  (2.8)  that 

VV*  =  -Y-  vR  .  (3.3) 

The  potential  component  4>p  is  a  homogeneous  solution  of  (2.8)  which  is  continuous  in  X 
and  satisfies  the  cascade  periodicity  and  acoustic  wave  propagation  requirements  stated 
above.  The  component  <j>R  is  a  particular  solution  of  (2.8)  (or  (3.2)  and  (3.3))  which 
satisfies  the  cascade  periodicity  requirement  and  the  wake  jump  conditions  (2.10)  and 
(2.11).  We  seek  first  a  solution  for  the  component  of  the  velocity  potential  describing 
the  acoustic  response  in  the  far  upstream  and  far  downstream  flow  regions  in  Section 
(3.2)  and  then  proceed  to  determine  a  solution  for  the  component  of  the  potential 
describing  the  response  to  vortical  waves  in  Section  (4.2).  The  potential  component 
4>R  depends  on  the  rotational  velocity  fluctuations  in  the  far  downstream  region,  and 
a  solution  describing  these  fluctuations  will  be  given  in  Section  (4.1).  We  will  use  the 
independent  variables,  (£,  tj)  and  (Tm,  Nm),  for  determining  the  required  solutions. 
The  Tm,  iVm-Cartesian  axes  are  parallel  and  normal,  respectively,  to  the  mth  wake  in 
the  far  downstream  region  and  have  their  origins  at  the  intersection  of  mth  wake  and 
the  line  £  =  £+  (see  figure  1).  Therefore  we  can  write 

Tm  +  iNm  =  [£  -  £+  +  i(r]  —  r?+  —  mG)\ exp(-t'fl0O)  .  (3.4) 

Unit  vectors  in  a  given  coordinate  direction,  say  the  ^-direction,  are  indicated  by  ef. 
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3.2  The  Potential  due  to  Acoustic  Response  Disturbances 

To  simplify  the  discussion  in  this  section  we  will  determine  4>P{X)  for  acoustic 
response  disturbances  only.  Once  this  solution  is  obtained,  it  becomes  a  simple  matter 
to  add  to  it  the  appropriate  terms  that  account  for  incident  acoustic  excitations.  Also, 
we  will  omit  the  subscripts  ^oo  when  referring  to  free-stream  flow  properties,  but  it 
must  be  recalled  that  the  symbols  V,  0  M,  etc.  represent  either  inlet  or  exit  free- 
stream  flow  properties.  Finally,  the  following  derivation  applies  to  subsonic  ( M  <  1), 
transonic  (in  the  limit  M  — »  1)  and  supersonic  (M  >  1)  flows.  However,  it  is  restricted 
to  flows  with  subsonic  axial  velocity  component,  i.e.,  M cos  fl  <  1.  Thus  if  M  >  1, 
|  tan  n|  >  (M2  —  l)1/2,  where  0  is  the  free-stream  flow  angle  and  ±tan  1[(M2  —  l)1/2] 
are  the  angles  between  the  flow  Mach  waves  (or  characteristic  lines)  and  the  free-stream 
flow  direction.  Once  this  derivation  is  completed,  we  will  simply  present  for  the  sake 
of  completeness  the  appropriate  results  for  supersonic  axial  flows  (M  cos  fi  >  1). 

The  desired  solution  for  the  unsteady  potential  component  <j)p  can  be  determined 
by  Fourier  methods.  Thus  we  consider  periodic  fundamental  solutions  of  the  form 

'finiZiV)  =  fn{  OeXPilkrimV),  n  =  0,  ±1,  ±2,  ...  ,  (3.5) 


where  we  use  the  symbol  k  to  denote  the  wave  number  of  a  response  fluctuation  and 
kn,n  =  (cr  +  2nn)G~1  is  a  wave  number  in  the  cascade  circumferential  or  redirection. 
After  substituting  (3.5)  into  the  far-fleld  differential  equation  (2.8)  and  solving  the 
resulting  ordinary  differential  equation  for  /„,  we  find  that 


MO 


“i.»  exp(xi,n£)  +  Oj.n  exp(xi,i.O  for  4/° 

(fll.n  +<l2,n£)eXp(Xl,n£)  f°r  4=° 


(3.6) 


where 


X  j  —  id^  T  tJVf  &n  cos  H  , 

2»n 


6n  =  (ojV  1  +  kn<n  sin  fi) /(l  —  M 2  cos2  H)  , 
d2n  =  (l-M2cos2n)-1A;2)n-M2«52 
and  dn  is  the  principal  root  of  d2  ;  i.e., 


f  |d„|,  if  d2  >  0 
(  i  jd„j,  if  d2  <  0 


(3.7) 


(3.8) 


If  d2  ^  0,  each  fundamental  solution  describes  two  wave-like  disturbances  which, 
depending  upon  the  sign  of  d2 ,  either  grow  or  decay  exponentially  in  the  axial  direction 
(d2  >  0)  or  propagate  carrying  energy  away  from  or  toward  the  blade  row  (d2  <  0).  The 
disturbance  which  shows  growth  with  increasing  axial  distance  from  the  blade  row  or 
propagation  of  energy  toward  the  blade  row  must  be  eliminated.  The  condition  d2  =  0 
which  divides  these  two  types  of  behavior  is  called  the  acoustic  resonance  or  “cut-off” 
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condition.  Here,  the  second  term  on  the  right-hand-side  of  (3.6)  must  be  eliminated 
because  the  unsteady  potential  must  remain  bounded  with  increasing  distance  from 
the  blade  row.  Thus  in  each  case,  d„  >  0,  dl<  0  or  d2n  =  0,  we  can  write 

(^i  *))  an  \{hndn  T  fiVf  6n  COS  n)  £  T  >  (3.9) 

and  our  task  is  to  determine  the  correct  value  of  hn  (i.e. ,  +1  or  —1)  for  d2n  ^  0. 

It  follows  from  (3.7)  and  some  algebra  that  the  circumferential  wave  numbers  at 
which  resonance  occurs,  k^  ,  are  given  by 

=  ojV~1M{  1  -  M2)~1{Msmn±\  1~M2  cos2  f I),  M  ^  1.  (3.10) 

Note  that  |Msinfi|  >  (1  —  M2  cos2  fl)1'2  for  M  >  1.  Also,  the  ranges  of  circumferential 
wave  numbers  for  which  disturbances  vary  exponentially  in  the  axial  direction  (i.e., 
d 2n  >  0)  are  given  by 


kn,n  < 

for  M  <  1 

and  uj  <  0  , 

(3.11a) 

kn,n  >  k ^ 

for  M  <  1 

and  uj  <  0 

(3.11b) 

and 

k <  k 

for  M  >  1 

and  u  <  0  . 

(3.12) 

Those  for  which  disturbances  propagate  in  the  axial  direction 

by 

k^  <  kn  n  <  k^  for  M  <  1  and  cu  <  0  , 

(i.e.,  d2n  <  0)  are  given 

(3.13) 

kn  n  >  k 

for  M  >  1 

and  uj  <  0 

(3.14a) 

and 

kn,n  <  k(n±] 

for  M  >  1 

and  u  <  0  . 

(3.14b) 

In  the  transonic  limit  ( M  l)  we  find  that 

[  2ujV~1M2  sin  f2/(l  —  M2)  +  u>V~l  /  {2  sin  fi)  ,  ft>0 
Jim  *{*>  -  (3.15) 

[  — cuC  1/(2sinfl)  ,  0  >  0 

and  the  results  achieved  using  the  present  analysis  will  agree  with  those  based  on  a  sim¬ 
ilar  analysis  of  the  high-frequency  transonic  small-disturbance  equation  (i.e.,  (2.8)  with 
Atoo  set  equal  to  V^).  The  high-frequency  transonic  equation  (Landahl  1961)  yields 
only  a  single  resonant  circumferential  wave  number,  i.e.,  k^=l^  =  — c<jV'~1/(2sin  H), 
and  disturbances  at  circumferential  wave  numbers  to  one  side  of  this  resonant  wave 
number  are  of  propagating  type  while  those  at  circumferential  wave  numbers  to  the 
other  side  are  of  decaying  type. 
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It  is  a  simple  matter  to  determine  the  appropriate  form  of  the  fundamental  solution 
for  tangential  wave  numbers  satisfying  the  inequalities  in  (3.11)  and  (3.12).  Since 
physically  admissible  solutions  must  be  bounded  in  the  far  field,  hn  must  be  positive 
in  the  far  upstream  region  (£  <  £_)  and  negative  in  the  far  downstream  region  (£  > 
£+).  For  tangential  wave  numbers  satisfying  the  inequalities  in  (3.13)  and  (3.14)  the 
fundamental  solution  (3.9)  describes  a  wave  that  propagates  in  the  axial-flow  direction 
at  wave  number 

k^n  =  hn  |  dn  |  +M26ncosfl  (3.16) 

and  the  determination  of  the  correct  value  of  hn  is  somewhat  more  difficult.  Various 
methods  have  been  proposed  to  dictate  this  choice.  These  involve  adding  an  artificial 
damping  term  (e^f)  to  the  left-hand-side  of  the  far-field  differential  equation  (2.8) 
(Verdon,  Adamczyk  &  Caspar  1975;  Adamczyk  1978),  considering  an  unsteady  excita¬ 
tion  with  slowly  growing  temporal  amplitude,  i.e.,  with  w  =  u)'  i A,  where  A  is  a  small 
positive  number  (Ni  1979),  or  determining  a  far-field  solution  in  terms  of  coordinates 
moving  at  the  uniform  inlet  or  exit  velocity  (Whitehead  &i  Davies  1983).  In  the  first  two 
methods  the  correct  form  of  the  fundamental  solution  is  determined  by  requiring  that 
the  unsteady  potential  remains  bounded  as  |£|  — *  oo.  Once  this  form  is  determined,  e<j 
or  A  are  allowed  to  approach  zero.  In  the  method  of  Whitehead  &:  Davies  (1983)  the 
form  of  the  fundamental  solution  in  the  moving  coordinate  frame  which  ensures  that 
acoustic  energy  propagates  away  from  the  blade  row  is  used  to  determine  the  correct 
value  of  hn  in  (3.9)  or  (3.16).  This  method  seems  to  have  the  most  physical  appeal  and 
is  therefore  the  one  adopted  herein. 

Consider  coordinate  axes  £'  and  t]'  parallel  to  the  space-fixed  £-  and  77-axes,  re¬ 
spectively,  but  moving  at  the  velocity  V .  After  replacing  <f>(X)  by  0(X,t)  =  4>(X)etui 
and  tix)  by  d/dt  in  the  far-field  equation  (2.8)  and  applying  the  independent  variable 
transformation 

£'  =  £  —  V$t,  r}'  =  rj  —  Vnt  and  t'  =  t  (3-17) 

we  find  that 

—  A  20t't>  =  0  .  (3.18) 

Equation  (3.18)  has  a  fundamental  solution  of  the  form 

4>{X,t')  =  F(X  ■  ep -At')  +  G(X  •  7P  +At')  ,  (3.19) 

where  X1  is  a  position  vector  in  the  £',  ^'-coordinate  frame  and  ep  is  a  unit  vector 
pointing  in  the  direction  of  wave  propagation.  Since  only  those  waves  which  propagate 
away  from  the  blade  row  are  admissible,  ep  •  e must  be  less  than  zero  for  £'  — *  —  oo 
and  greater  than  zero  for  £'  — >■  oo.  Then,  since  the  function  G  describes  waves  that 
propagate  toward  the  blade  row,  it  must  be  set  equal  to  zero,  and  the  nth  term  of  the 
general  solution  to  (3.18)  can  be  written  as 

=  b'nexp[i(x)'n(t'  -  £'cos0n/A  -  rj' sin  d„/ A]  .  (3.20) 
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This  solution  describes  a  wave  propagating  at  speed  A  and  at  angle  0n,  measured 
counterclockwise  from  the  £'-axis,  where  7r / 2  <  6n  <  37t/2  for  £'  — *  —  oc  and  \6n\  <  7r/2 
for  £'  — *  oo.  The  temporal  frequency  of  the  nth  disturbance  wave  as  seen  by  an  observer 
moving  with  the  mean  fluid  velocity  V  is  u)'n. 

A  comparison  of  the  result  in  (3.20)  with  the  fundamental  solution  for  d?n  <  0  in 
(3.9)  dictates  the  correct  value  of  hn  to  be  used  in  the  latter  equation.  It  follows  from 
(3.7),  (3.9),  (3.16),  and  (3.20)  that 

u'n  =  w+  kn  ■  V=  V(6n  +  hn  |  dn  |  cos  fi)  , 

Sin 

and 

cos6n  =  -kt  n  ,  (3.21) 

where  kn>n  satisfies  one  of  the  inequalities  in  (3.13)  and  (3.14).  Now,  a  disturbance 
carried  by  the  nth  wave  travels  in  the  axial  direction  with  velocity  C^n,  where 

C^n  —  Vj  +  A  cos  0n  =  —hnA2(  1  —  JVf2  cos2  U)  \  dn  \  /u>'n  .  (3.22) 

The  right-hand-side  of  (3.22)  follows  from  (3.16)  and  (3.21).  Note  that  in  the  limit 
as  d2n  — >  0  (i.e.,  as  a  resonance  condition  is  approached),  C^n  — ♦  0  and  there  is  no 
transfer  of  acoustic  energy  in  the  axial  direction.  Since  the  velocity  Cc  n  must  be  less 
than  zero  in  the  far  upstream  region  and  greater  than  zero  in  the  far  downstream 
region,  the  correct  value  of  hn  in  (3.16),  (3.21)  and  (3.22)  depends  on  the  sign  of  oj'n. 
Some  analysis  of  (3.7),  (3.8),  (3.13)  and  (3.14)  shows  that  Sn  <  ±|dn|cosf2  for  oj  <  0 
if  the  free-stream  flow  is  subsonic  and  k^'  <  kn  n  <  k^\  or  if  the  free-stream  flow  is 
supersonic  with  the  product  wfl  >  0  and  kn  n  >  k^  or  with  cuff  <  0  and  k^  n  <  kfl 
Under  these  conditions  c ofn  and  to  have  the  same  sign;  therefore,  hn  =  sgfr?(cj)  in  the  far 
upstream  (£<£_)  region  and  hn  =  —  sgn(uj)  in  the  far  downstream  (f  >  f+)  region, 
where  sgn(w)  —  ±1  for  u  <  0.  If  the  free-stream  flow  is  supersonic  with  uU  >  0  and 
kn,n  <  or  with  toft  <  0  and  n  >  k^\  then  6n  >  ^\dn\  cos  f]  for  uo  <  0.  In  this  case 
uj'n  and  to  are  of  opposite  sign  and  hn  —  =Fsgn(u)  for  £  >  . 

Thus,  as  a  result  of  the  foregoing  analysis,  we  can  write  the  solution  for  the  unsteady 
potential  in  the  far  field  due  to  the  acoustic  response  fluctuations  as 

_ >  oo 

4>P{X)  =  Y,  a",T°°  exP[Xn(£  -  t?)  +  ikrtA'n  -  »7t)],  t  >  Zt  ■  (3.23) 

n  —  ~  oo 

The  constants  an i=F00  are  determined  by  matching  these  analytic  solutions  to  a  near-field 
(numerical)  solution,  (^_,r?_)  and  (£+,r/+)  are  reference  locations,  i.e.,  the  points  of 
intersection  of  the  reference  blade  (m  =  0  )  stagnation  streamline  and  the  lines  f  =  £_ 
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and  £  =  £+,  respectively,  and 


f  ±  1 

dn 

\+iM26n  cos  n  I 

(3.24a) 

Xn  =  Pn+ike,n  =  \  *(  ±Sgn(w) 

K 

|+M2<5„cosfl)  > 

for  £  >  ,  (3.24b) 

{  *(  Tsgn(w) 

I  dn 

j  +M26n  cos  f2)  J 

(3.24c) 

where  dn  and  S„  are  defined  in  (3.7)  and  (3.8).  At  resonance  (i.e.,  knn  —  fcW  (see 
(3.10)))  |d„|  =  0.  If  \dn\  0,  Xn  is  determined  using  (3.24a)  if  the  inequalities  in 
(3.11a),  (3.11b)  or  (3.12)  are  satisfied;  (3.24b)  if  the  inequalities  in  (3.13),  (3.14a)  for 
cufi  >  0  or  (3.14b)  for  u>U  <  0  are  satisfied;  and  (3.24c)  if  the  inequalities  in  (3.14b) 
for  u)U  >  0  or  in  (3.14a)  for  wfl  <  0  are  satisfied.  This  information  is  summarized  in 
figures  2  and  3.  In  the  limit  M  — >  1  the  subsonic  (see  figure  2)  and  supersonic  (figure  3) 
results  become  identical,  and  therefore  no  disturbances  of  the  type  3.24c  can  exist  in  a 
transonic  free  stream.  It  is  also  important  to  note  that  the  expression  for  the  potential 
due  to  prescribed  far  upstream  and/or  far  downstream  acoustic  excitations  must  have 
the  same  form  as  (3.23)  except  that  the  Xn  must  be  replaced  by  —  xn- 

An  analysis  similar  to  the  above  can  be  performed  for  supersonic  relative  flow  with 
supersonic  axial  velocity  component  (see  also  Ni  1979;  Whitehead  1987).  In  this  case 
the  fundamental  solutions  for  4>p  are  also  given  by  (3.5),  (3.6)  and  (3.7),  but  d2n  is  less 
than  zero  for  all  circumferential  wave  numbers.  Therefore  each  fundamental  solution 
describes  two  propagating  acoustic  response  disturbances.  It  is  easy  to  show  that 
both  of  these  disturbances  carry  energy  in  the  downstream  axial  direction.  Hence, 
for  supersonic  streams  with  M  cos  fl  >  1  the  potential  0P  due  to  acoustic  response 
disturbances  has  the  form 


<j>p(X)  = 


o  ,  £  <  £- 

oo 

2  {Gi,nexp[xi,„(£  -  £+)]  +  a2,nexp[x2.n(£  -  £+)]} 

=  —  oo 

x  exp[ik„<n(r]  -  ??+)]  ,  £  >  £+  , 


(3.25a) 

(3.25b) 


where  Xi.n  and  X2,n  are  defined  by  (3.7)  and  (3.8). 

In  a  supersonic  stream  with  supersonic  axial  velocity  component  only  far  upstream 
acoustic  excitations  can  affect  the  flow  through  the  blade  row.  The  potential  due  to  such 
excitations  has  the  same  form  as  (3.25b),  except  that  £+  and  r/+  must  be  replaced  by  £_ 
and  ri_ ,  respectively.  Because  of  the  manner  in  which  unsteady  disturbances  propagate 
in  supersonic  streams,  an  analytic  solution  for  the  velocity  potential  response  in  the 
far  upstream  or  far  downstream  region  is  not  required  to  determine  the  unsteady  flow 
through  a  cascade  operating  at  supersonic  inlet  or  exit  velocity.  Thus,  our  concern 
here  is  primarily  with  cascades  operating  at  subsonic  axial  inlet  and  exit  velocities  for 
which  a  knowledge  of  the  far-field  potential  response  is  indeed  required  to  predict  near 
field  unsteady  flows.  Hence,  throughout  the  remainder  of  this  report  when  we  refer 
to  supersonic  inlet  or  exit  conditions,  it  is  assumed  implicitly  that  the  axial  velocity 
component  is  subsonic. 
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4.  THE  UNSTEADY  POTENTIAL  DUE  TO  VORTICAL 

DISTURBANCES 


The  unsteady  potential  component  <f>R  is  determined  as  a  solution  of  the  differential 
equations  (3.2)  and  (3.3)  which  satisfies  cascade  periodicity  (2.4)  and  the  conditions  of 
continuity  of  pressure  (2.10)  and  normal  velocity  component  (2.11)  across  blade  wakes. 
In  one  sense  it  is  a  much  simpler  task  to  determine  <j>R  than  it  is  to  determine  the 
far-field  potential  (f)p  because  a  single  expression  holds  for  (j)R  regardless  of  the  sign 
of  M  —  1,  M  cos  f7  —  1,  fi  or  u.  However,  the  potential  component  <f>R  depends  upon 
the  rotational  velocity  fluctuation  and  therefore,  we  must  first  determine  a  solution  for 
vff(£  >  £+,f?)  before  proceeding  to  determine  the  corresponding  solution  for  (f>R . 

4.1  The  Rotational  Velocity  in  the  Far  Downstream  Region 

The  rotational  velocity  in  the  far  downstream  region  is  determined  as  a  solution  to 
the  field  equation  (2.9)  for  a  periodic  (in  r?)  distribution  of  rotational  velocity  along  £  — 
£+  which  can,  in  general,  be  discontinuous  at  blade  wakes.  To  determine  this  solution 
it  is  again  convenient  to  decompose  the  dependent  variable  into  two  components.  Thus 
we  set 

+  Vd  >  l4*1) 

where  vR  is  continuous  and  vR  possesses  jump  discontinuities  at  blade  wakes,  i.e., 

-  uV^Tm)}  on  Wm,  m  =  0, ±1, ±2, ...  .  (4.2) 

We  require  that  both  vR  and  vR  satisfy  (2.9)  and  the  periodicity  condition  (2.4)  and 
that  V  *  Vd  “  0-  It  then  follows  from  (2.9)  and  (2.4)  that 

R  00 

vc=  bn  exp[i  kn  -(X  -X+)],  £  >  £+  ,  (4.3) 

ft  =  —  oo 

where  the  Fourier  coefficients,  6n,  are  determined  in  terms  of  the  rotational  velocity 
distribution  on  f  (see  (4.24)  below)  and 

hn  —  (t^oo  SeC  floo  T  tan  floo)^ 

=  -wV^'er  +  (wF'Hanfioo  +  kn.n secU^e N  ,  (4.4) 

where  kn  n  =  (a  +  2nn)G~1. 

A  particular  solution  for  the  discontinuous  component  of  the  rotational  velocity, 
v$ ,  in  the  far  downstream  region  of  the  mth  passage  (i.e.,  in  the  region  Tm  >  0,  0  < 
Nm  <  Gcosfioo)  can  be  determined  conveniently  in  terms  of  the  Cartesian  coordinates 
Tm  and  Nm.  Once  this  result  is  obtained,  it  can  be  extended  to  provide  a  solution  in 
terms  of  the  independent  variables  £  and  which  holds  over  the  entire  downstream 
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region.  It  follows  from  (2.9)  and  (4.2)  that  the  required  solution  for  the  mth  passage 
can  be  written  in  the  form 

vRd=Fm  (iVm)exp(-iwF~]rm)  .  (4.5) 

Since  vd  is  divergence  free  and  has  a  jump  discontinuity  across  blade  wakes, 

-fwV-1  Fm-7T  +  Fm  ■  7N=  0  (4.6) 

and 

Fm  (0+)-  Fm  (Geos  0”)  exp[— i(c  -f  wV^Gsin  0)]  =  |v  Jfle/exp(fma)  ,  (4.7) 

respectively.  A  solution  for  Fm  which  meets  the  requirements  (4.6)  and  (4.7)  can  be 
constructed  as  a  linear  combination  of  the  exponential  functions  exp(u)Vr(x;1Arm)  and 
exp(— uV^Nm)-  After  making  several  trials  and  performing  some  algebra,  we  find 
that 

Fm  ( Nm )  =  ew[C1exp(u;I^1iVm)(er  -f*  eN) 

+  G2 exp(— wV;0^1iVm)( ex  -t  eN)}  ,  (4.8) 

where  the  constants  C\  and  G2  are  given  by 

Gj  =  hi  -  exp(wV0^1Ge',n°°  -  ia)]|v  ]i?e/  •  (er  -t'cjv)  (4.9a) 

and 

G2  =  -[1  -  exp(-wFo;1Ge,no0  -  ia)lvRjFef  •  (‘ eT  +^n)  ■  (4.9b) 

2 

The  solution  for  the  rotational  velocity  over  the  entire  far  downstream  region  (£  > 
£+)  is  then  determined  by  combining  the  results  in  (4.1),  (4.3),  (4.5),  (4.8)  and  (4.9) 
to  obtain 

vR  (X)  =  jr  bn  exp[t  -(X  -  X+)] 

n=  — oo 

+  C\  exp[— (X  —  X_{-)*(e£  +  ter7)e  00  ]  ( e  £  e  n)e 

oo 

x  exp[imo  -  Tnu>V~1Ge~'Qcc}[U(ri  -  riWm)  -  u(r)  -  r)Wm+l)} 

m=— oo 

+  G2exp[-twF-1(X  -  X+)  •  (e€  -i  e„)e,n~](e<  -i  e„)e*n~ 

OO 

x  Y  exp[un<r  +  rrvJjV-lGem~}\U{ri  -  rjWm)  -  U(t]  -  vWm+l)}  •  (4-10) 
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Here 

V  =  VWm  =  h+  +  mG  +  (£  -  £+)  tan  (4.11) 

defines  the  location  of  the  straight-line  approximation  to  the  mth  wake, 


x-x+=  ({-{+)?( +(,->,+)  7, 

=  ( Tm  +  mG  sin  er  +{Nm  +  mG  cos  O^)  "ejv 


(4.12) 


and 


e  j  ±1  (~eT  ±i  ~eN)  exp(±?noo)  .  (4.13) 

The  divergence  of  the  rotational  velocity  in  the  far  downstream  region  is  given  by 


V-  v  =  ^  i  kn  ■  bn  exp[t  kn  -(X  - 


(4.14). 


4.2  The  Potential  due  to  Vortical  Disturbances 

As  in  the  case  of  the  rotational  velocity,  we  seek  first  a  solution  for  the  potential 
component  <j>R  in  the  far  downstream  region  of  the  mth  passage  (i.e.,  Tm  >  0,0  < 
Nm  <  G  cosfl)  and  then  extend  this  result  to  provide  a  solution  over  the  entire  far 
downstream  region  (f  >  f+).  It  follows  from  the  field  equations  (3.2)  and  (3.3),  the 
cascade  periodicity  (2.4)  and  wake-jump  conditions  (2.10)  and  (2.11),  and  the  foregoing 
results  for  the  far  downstream  rotational  velocity  that  the  required  solution  must  be 
of  the  form 

<f>R(Tm,  Nm)  =  ifm(ATm)exp(-tcuV“1rm)  ,  (4.15) 

where 

OO 

Hm{Nm)  -  u>2V^2H(Nm)  =  -  J2  i  kn  ■  t>n  exp {i[kN>nNm  -f  m(cr  +  27m)]}  ,  (4.16) 

n=  — oo 

Hm{ 0  )  —  exp[— t(a  +  wV^^sin  noo)]//(Gcos  n^,)  =  exp(tmCT)  (4.17) 

and 

H'm{0+)-exp\-i((7+ujVoolG smCloo^H^G cosfl^)  =  —  \v  e^exp(fm(r).  (4.18) 

The  general  solution  of  (4.16)  is 

Hm[Nm)  =  Zh  exp(tma)  exp^V^1  Nm)  +  D2  exp(tma)  exp(-c vV~lNm) 


^  i  k  ■  b 

+  ^ - -exp{i[kN>nNm  + m(a +  2nn)]}  .  (4.19) 

n=-°°  |^n|2 
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The  constants  D\  and  Z?2  are  evaluated  using  the  wake  conditions  (4.17)  and  (4.18) 
and  are  given  by 

Di  =  1[1  -  exp(— 1<7  +  wV-1Ge_,'noo)]"1(|^]li?e/  -  «w)  (4.20a) 

and 

D2  =  2[l-exp(-ta-wFro1Ge,'n°o)j"1(|(/)l/?e;  +  w~1V00[v/?Jfie/-  eN)  .  (4.20b) 

Therefore,  the  solution  for  the  velocity  potential  fluctuations  in  the  far  downstream 
region  of  the  mth  passage  due  to  rotational  velocity  disturbances  and  wake  vorticity, 
is  given  by 


<PR{Tm,  Nm)  =  ei””{Dhexp[-i«,V-1{Tm  +  iNm)\  +  D2exp[-iuV~1{Tm  -  *JVm)]} 
i  k  ■ b 

+  E  — =? — -  exp{i[kTjnTm  +  kN}nNm  +  m((j  +  2mr)}} 
n=-°°  |  kn  |2 

Finally,  the  solution  for  <f)R  over  the  entire  downstream  region  £  >  £+  is 
4R(X)  =  Dr  exp [-iuV-'iX  -  X+)  ■  (e€  +i  en)e-n-] 

OO 

x  E  exp\ima-u)V~imGe-iU^}[U{r]-ilwm)-U{r)-r]Ww+l)} 

771—  —  OO 

+  D2exp\-iuV~1(X  -X+ )  •  (e*  -i  e„)e''n~] 


(4.21) 


x  E  exp[tma  +  uj'VjmGe^lUit]  -  i]Wtn)  ~U(r)  -r?B.m+J)] 

771—  —  OO 

°°  7  k  .  7  _ 

+  E  "  exP[*  kn  \X  -  X+)}  .  (4.22) 

n=-°°  |^n|2 

4.3  The  Complete  Solution  for  the  Unsteady  Potential 

We  have  considered  the  flow  through  an  isolated  two-dimensional  blade  row  in  which 
unsteady  fluctuations  are  produced  by  small-amplitude  periodic  (in  77  and  t)  structural 
(blade  motions)  and/or  external  aerodynamic  excitations  (incident  entropic,  vortical 
or  acoustic  disturbances).  Far  from  the  blade  row  the  complex  amplitude,  <f)(X),  of  the 
linearized  unsteady  potential  can  be  expressed  as  the  sum  of  two  components  (3.1) — 
one,  <j>p(X),  associated  with  acoustic  disturbances  and  the  other,  <j>R( X),  associated 
with  vortical  disturbances.  Analytic  solutions  for  <j)P  and  <pF  due  to  acoustic  and 
vortical  response  fluctuations  are  given  in  (3.23)  for  subsonic  axial  flow  and  (4.22), 
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respectively.  The  unsteady  potential  component  <f>R  depends  upon  the  far  downstream 
rotational  velocity  field  which  is  described  in  (4.10). 

The  physical  variables  that  the  potential  component  (f)P  depends  on  explicitly  are 
the  free-stream  Mach  number  and  flow  angle,  the  cascade  blade  spacing,  the  frequency 
and  interblade  phase  angle  of  the  unsteady  excitation  and  the  complex  amplitudes  of 
the  first-harmonic  acoustic  disturbances.  In  particular,  the  complex  amplitudes,  an)T0 0, 
of  the  acoustic  response  fluctuations  are  given  by 

rri^+G 

fln.Too  =  G  /  [<£(£-,77)  -  4>R{Z+,v)U(£  ~  C+)}exp(-ikninrj)dji  .  (4.23) 

Note  that  if  a  trailing-edge  Kutta.  condition  is  imposed  to  determine  the  steady  flow 
through  the  blade  row,  the  exit  free-stream  conditions  cannot  be  prescribed  arbitrarily, 
but  are  determined  by  the  conditions  imposed  at  the  cascade  inlet.  The  physical 
variables  that  the  potential  component  (f>R  depends  on  are  the  far  downstream  free- 
stream  flow  speed  and  angle,  the  blade  spacing,  the  frequency  and  interblade  phase 
angle  of  the  unsteady  excitation,  the  complex  amplitudes  of  the  rotational  velocity 
fluctuations,  6n,  where 


— *  rri++G  

bn=G~  /  [v  (£+>*?)-  vd  (£+,?7)]exp(-^n77)d?/  ,  (4.24) 

and  the  jumps  in  the  potential  and  rotational  velocity  at  the  intersection 

of  the  reference  wake  and  the  axial  line  f  =  £+.  Note  that  if  the  rotational  velocity 
is  continuous  across  blade  wakes,  the  constants  C\  and  C2  in  (4.10)  are  identically 
zero  and  the  constants  D\  and  D2  in  (4.22)  depend  only  on  the  jump  in  potential 
M Ref-  In  addition,  if  there  are  no  prescribed  incident  rotational  velocity  disturbances, 
vR{X)  ~  0  throughout  the  field  and  the  last  term  in  (4.22)  is  identically  zero.  In  this 
important  special  case  (f>R  is  driven  only  by  the  vorticity  shed  at  the  blade  trailing 
edges  and  convected  along  the  blade  wakes,  which  is  directly  related  to  the  fluctuating 
aerodynamic  loads  acting  on  the  blades. 

4.4  Implementation  of  Far-Field  Solutions  in  Numerical  Computations 
The  complex  amplitudes  an,^oc  and  the  potential  and  rotational  velocity  jumps 
and  | and  hence,  the  far-field  solutions  for  <f>  and  vR  are  known  in  terms  of 
the  potential  distributions  along  the  lines  £  =  and  the  rotational  velocity  distribu¬ 
tion  along  £  ~  £+.  This  feature  can  be  applied  to  construct  a  discrete  approximation  to 
the  linearized  unsteady  boundary-value  problem  over  an  extended  blade  passage  region 
bounded  by  the  axial  lines  £  =  .  For  example,  if  the  computational  mesh  consists  of 

one  set  of  lines  parallel  to  the  blade  row  (see  Verdon  &  Caspar  1982,1984),  then  the 
discrete  approximation  to  the  field  equation  (2.8)  at  mesh  points  on  the  boundary  lines 
£  =  £=f  depend  upon  values  of  <fr  and  vR  upstream  and  downstream,  respectively, 
of  these  axial  mesh  lines.  However,  because  the  potential  upstream  of  £  ~  £  and  the 
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potential  and  rotational  velocity  downstream  of  £  =  £+  are  known  (see  (3.23),  (4.10) 
and  (4.22))  in  terms  of  the  potential  distribution  on  £  =  £_  and  the  potential  and 
rotational  velocity  distributions  on  £  =  £+,  the  discrete  approximations  to  the  field 
equations  on  £  =  can  be  constructed  entirely  in  terms  of  0  and  vR  values  at  points 
on  the  upstream  and  downstream  boundaries  and  in  the  interior  (£_  <  £  <  £+)  of  the 
computational  domain.  Thus,  using  the  analytic  far-field  solutions,  we  can  write  a  dis¬ 
crete  equation  for  each  unknown  <p  value  which  depends  only  upon  information  within 
and  on  the  boundaries  of  the  computational  domain.  Hence,  the  system  of  discrete 
equations  is  completely  specified. 

The  foregoing  approach  has  been  used  in  the  blade  flutter  calculations  (s  —  vR  = 
pj  =  0)  (see  Whitehead  1982;  Verdon  &;  Caspar  1982,  1984;  Usab  &;  Verdon  1986)  and 
is  a  particularly  convenient  one  for  cascades  operating  under  subsonic  inlet  and  exit 
conditions  because  only  a  limited  number  of  acoustic  response  waves  persist  in  the  far 
field.  For  supersonic  flows  with  subsonic  axial  velocity  component,  the  Fourier  series 
representation  for  <j)P  (3.23)  may  not  be  so  useful  because,  in  principle,  many  terms  of 
this  series  will  be  needed  to  represent  the  abrupt  changes  in  the  velocity  potential  that 
occur  at  Mach  waves.  We  should  also  note  that,  in  addition  to  their  use  in  constructing 
complete  sets  of  finite-difference  equations,  the  understanding  of  the  unsteady  flow  in 
the  far  field  gained  through  an  examination  of  the  analytic  far-field  solutions  is  essential 
to  an  accurate  resolution  of  the  flow  through  the  cascade.  In  particular,  mesh  spacings 
must  be  selected  so  that  the  important  variations  in  the  rotational  velocity  and  the 
velocity  potential,  indicated  by  the  unsteady  flow  behavior  in  the  far  field,  can  be 
resolved  accurately. 
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5.  UNSTEADY  FLOW  BEHAVIOR  IN  THE  FAR  FIELD 


We  proceed  to  analyze  the  unsteady  aerodynamic  response  in  the  far  field  of  an 
isolated  two-dimensional  cascade  operating  at  subsonic  axial  inlet  and  exit  velocity  by 
examining  the  parametric  dependence  of  the  wave  numbers  and  attenuation  constants, 
that  control  the  spatial  variation  of  far-field  potential  disturbances.  This  type  of  infor¬ 
mation  is  a  prerequisite  for  prescribing  the  computational  mesh  needed  for  a  numerical 
resolution  of  the  unsteady  flow  through  the  blade  row.  Thus  we  will  examine  the  be¬ 
havior  of  the  velocity  potential  in  the  far  upstream  and  downstream  regions  caused  by 
acoustic  response  disturbances  and  the  behavior  of  the  potential  in  the  far  downstream 
region  caused  by  continuous  rotational  velocity  disturbances  and  concentrated  vorticity 
convected  along  the  blade  wakes.  But  first  we  will  provide  a  useful  classification  for 
unsteady  excitations  which  is  based  upon  the  acoustic  responses  that  they  produce  in 
the  far  field. 

5.1  Classification  of  Unsteady  Excitations 

We  consider  a  fundamental  pressure  response  disturbance  and  using  (3.10)  define 
the  constant  Cn  as 

_  \u>\M{  1  -  M2  cos2  0)s 

"  |(1  -M2)Vk^n  -  toM2smd\  1  '  ' 

where  the  subscripts  — oo  or  +oc  must  be  applied  to  M,  V  and  f 1  and  therefore  Cn  to 
refer  to  the  far  upstream  or  far  downstream  flow  regions.  We  then  classify  the  nth  pres¬ 
sure  response  disturbance  as  subresonant,  resonant  or  superresonant  depending  upon 
whether  Cn  is  less  than,  equal  to,  or  greater  than  one,  respectively  (see  also  Samoilovich 
1967;  Verdon  1979;  Bendiksen  1986).  For  a  prescribed  Mach  number,  flow  angle  and 
excitation  frequency,  resonance  ( Cn  —  1)  occurs  at  two  different  circumferential  wave 
numbers,  i.e.,  k^\  The  nth  acoustic  response  disturbance  is  superresonant  ( Cn  >  1) 
for  circumferential  wave  numbers  lying  between  these  two  resonant  values  (see  (3.12) 
and  (3.13))  and  subresonant  ( Gn  <  1)  for  circumferential  wave  numbers  less  than  the 
lower  and  greater  than  the  higher  of  the  two  resonant  circumferential  wave  numbers 
(see  (3.11)  and  (3.14)).  This  information  is  also  illustrated  in  figures  2  and  3. 

For  subsonic  (M  <  1)  relative  flow  k and  k M  have  different  signs.  Thus  one  res¬ 
onance  condition  (i.e.,  k^  for  u  <  0)  corresponds  to  a  response  disturbance  travelling 
in  the  positive  ^-direction  and  the  other  (i.e.,  k W  for  u  <  0)  to  a  wave  traveling  in 
the  negative  //-direction.  In  a  subsonic  free  stream  subresonant  disturbances  attenuate 
with  increasing  distance  from  the  blade  row  (i.e.,  Xn  is  determined  using  (3.24a))  while 
superresonant  disturbances  propagate  away  from  the  blade  row  (yn  is  determined  using 
(3.24b)).  It  follows  from  (5.1)  that  Cn  can  be  greater  than  one  only  for  a  limited  num¬ 
ber  of  acoustic  response  fluctuations  and  hence,  in  subsonic  flow,  that  such  fluctuations 
must  be  either  all  or  mostly  all  of  decaying  type. 

For  supersonic  flow  with  subsonic  axial  velocity  component  both  resonant  circum¬ 
ferential  wave  numbers  describe  waves  traveling  in  the  same  direction,  i.e.,  the  negative 
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r/'direction  if  both  u  and  fi  have  the  same  sign  and  the  positive  77-direction  if  u  and  0 
differ  in  sign.  In  a  supersonic  free-stream  superresonant  disturbances  attenuate  with 
increasing  axial  distance  from  the  blade  row  (i.e.,  Xn  is  determined  using  (3.24a)),  while 
subresonant  disturbances  propagate  and  carry  energy  away  from  the  blade  row  (xn  is 
determined  using  (3.24b)  or  (3.24c)).  Since  Cn  can  be  greater  than  one  only  for  a 
limited  number  of  acoustic  response  disturbances,  the  latter  are  either  all  or  mostly  all 
of  propagating  type. 

We  can  classify  unsteady  excitations  based  on  the  acoustic  response  that  they  pro¬ 
duce  in  the  far  field  as  follows.  Recall  that  the  far  upstream  and  far  downstream  free- 
stream  conditions  will  generally  differ.  Indeed,  the  flow  may  be  subsonic  in  one  far- field 
region  and  supersonic  in  the  other.  We  refer  to  an  unsteady  excitation  as  subresonant 
if  it  produces  only  acoustic  response  disturbances  which  attenuate  with  increasing  dis¬ 
tance  from  the  blade  row  at  subsonic  free-stream  Mach  numbers  or  only  disturbances 
which  persist  and  carry  energy  away  from  the  blade  row  at  supersonic  free-stream  Mach 
numbers.  A  resonant  excitation  is  one  for  which  at  least  one  response  wave  persists  and 
carries  energy  parallel  to  the  blade  row  in  either  the  far  upstream  or  far  downstream 
region  regardless  of  whether  the  free-stream  Mach  numbers  are  subsonic  or  supersonic. 
Finally,  a  superresonant  (m,n)  excitation  is  one  for  which  m,n  waves  persist  and  carry 
energy  away  from  the  blade  row  at  subsonic  inlet,  exit  Mach  numbers  or  attenuate  with 
increasing  distance  from  the  blade  row  at  supersonic  inlet,  exit  Mach  numbers.  Note 
that  if  either  m  or  n,  but  not  both,  is  zero,  the  excitation  is  regarded  as  superresonant. 
Thus,  if  the  inlet  and  exit  Mach  numbers  are  both  subsonic,  an  excitation  classified  as 
superresonant  (1,0)  produces  one  wave  which  carries  energy  away  from  the  blade  row 
in  the  upstream  axial  direction,  but  all  other  acoustic  response  disturbances  associated 
with  this  excitation  attenuate  with  increasing  axial  distance  from  the  blade  row. 


5.2  The  Acoustic  Response  in  the  Far  Field 

To  illustrate  the  acoustic  response  in  the  far  field  in  more  detail  we  introduce  the 
scaled  variables 


tn  =  u-'Vk^n  -  Vw~lG~l{o  +  27Tn) 

,  (5.2) 

=  u~lVIm{xn} 

(5.3) 

and 

P  =  oj-'V  \  Re{Xn}  | 

(5.4) 

where  Xn  is  defined  in  (3.24),  n  =  0, 1, 2, . . .,  and  we  assume 

that  the  interblade  phase 

angle  of  the  unsteady  excitation  varies  between  —n  and  +7 r,  i.e.,  ae[— 7r,7r’j.  The  res¬ 
onant  values  (i.e.,  t^)  of  the  scaled  circumferential  wave  number  depend  upon  the 
free-stream  Mach  number  and  flow  angle  (see  3.10),  while  the  axial  wave  number, 
£j,  and  the  attenuation  constant,  (3 ,  can  be  expressed  as  functions  of  the  free-stream 
Mach  number,  the  free-stream  flow  angle  and  the  tangential  wave  number  (see  (3.24) 
and  (3.7)).  Results  illustrating  the  behavior  of  these  parameters  for  selected  subsonic 
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(M  <  1)  flows  and  supersonic  ( M  >  1)  flows  with  subsonic  axial  velocity  compo¬ 
nent  (Mcosfi  <  1)  are  presented  below.  Although  these  results  depict  If.  and  0  as 
continuous  functions  of  it  must  be  emphasized  that  a  prescribed  excitation  at  tem¬ 
poral  frequency  u>  produces  only  a  discrete  set  of  acoustic  response  disturbances  at 

ln  =  (wG)_1F(< 7  +  27m), n  =  0,1,2, _  Also,  although  upstream  and  downstream 

acoustic  response  properties  are  depicted  below  for  the  same  free-stream  Mach  number 
and  flow  angle,  the  far  upstream  and  far  downstream  free-stream  states  will  generally 
differ. 

The  subsonic  results  and  the  supersonic  resonance  curves  presented  below  hold  for 
both  positive  and  negative  excitation  frequencies  (w  <  0)  and  non-negative  flow  angles 
(fi  >0);  the  supersonic  results  for  £^  and  0  hold  for  uj  >  0  and  fi  >  0.  Corresponding 
results  for  other  combinations  of  oj  and  fi  can  be  readily  determined  from  those  provided 
here.  Thus,  resonance  curves  for  17  <  0  are  obtained  by  changing  fi,  £j,  1  and  tn 
in  figures  4  and  7  to  -fi,  £ W  and  -Cn,  respectively.  Predictions  of  and  0  for 
M  <  1  and  fi  <  0  and  for  M  >  l,w  >  0  and  fi  <  0  can  be  obtained  from  those  given 
in  figures  5,  6,  8  and  9  by  changing  £,,  to  —  £„.  Predictions  of  £^  and  0  for  M  >  l,w  <  0 
and  fi  >  0  are  obtained  by  changing  i^-oo  and  £^i0o  to  C^)00  and  £|t  ,-oo>  respectively, 
in  figures  8  and  9.  Finally  predictions  of  £{  and  0  for  M  >  l,w  <  0  and  fi  <  0  are 
obtained  by  changing  Ln,  £^_oo  and  £^i00  in  figures  8  and  9  to  — £,,,  £^i0O  and  £j,-oo5 
respectively. 

5.2.1  Subsonic  Flow 

The  resonant  values  of  are  plotted  versus  Mach  number  (as  ordinate)  for  various 
free-stream  flow  angles  in  figure  4.  The  resonant  circumferential  wave  number  £|)1"-  varies 
substantially  with  Mach  number  and  flow  angle,  particularly  at  high  subsonic  Mach 
number,  but,  except  at  low  flow  angles,  the  dependence  of  £^  on  M  and  fi  is  rather 
limited.  Indeed,  as  M  — *  1  from  below,  £|j+^  — >  00  and  £ 1  — >  — (2sinfi)_1.  As  indicated 
in  figure  4,  the  extent  of  the  superresonant  region  or  interval,  i.e.,  (£|J"), £(,+^),  increases 
with  increasing  free-stream  Mach  number  and  to  a  lesser  degree  with  increasing  flow 
angle. 

For  example,  consider  an  unsteady  flow  driven  by  a  unit  frequency  excitation  (w  = 
1)  through  a  cascade  with  unit  blade  spacing  ( G  =  1)  and  operating  in  a  free  stream 
with  V  =  1  and  fi  =  45  deg.  These  conditions  are  somewhat  representative  of  those 
at  which  jet  engine  fan  or  compressor  blades  have  experienced  subsonic  flutter.  The 
calculations  leading  to  the  curve  for  fi  =  45  deg  in  figure  4  indicate  that  for  a  free- 
stream  Mach  number  of  0.7  and  interblade  phase  angles  lying  in  the  range  [ — 7r,  7tI 
all  acoustic  waves  attenuate  in  the  far  field  if  the  interblade  phase  angle,  er,  of  the 
unsteady  excitation  satisfies  the  conditions  a  <  —0.513  or  o  >  1.872,  but  that  one  such 
wave  (i.e.,  the  n  =  0  wave)  will  persist  and  carry  energy  away  from  the  blade  row  if 
—0.513  <  a  <  1 .872.  At  a  Mach  number  of  0.8  the  corresponding  ranges  are  a  <  —0.575 
or  a  >  3.09  and  0.575  <  a  <  3.09,  respectively.  The  situation  becomes  more  interesting 
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at  M  —  0.9  where  propagating  waves  occur  for  all  values  of  a.  In  particular,  these  waves 
occur  at  circumferential  wave  numbers  kntn  =  a  +  27rn  lying  in  the  range  (—0.639, 
6.668).  Both  the  n  =  0  and  n  =  1  acoustic  response  disturbances  persist  in  the  far 
field  if  —0.639  <  a  <  0.385,  but  only  the  n  =  1  or  the  n  =  0  disturbance  persists  if 
—  7r  <  a  <  —0.639  or  0.385  <  o  <  7r,  respectively.  As  the  Mach  number  is  increased 
further  for  our  example  configuration,  additional  acoustic  response  disturbances  will 
be  of  propagating  type.  Note  that  in  addition  to  the  dependence  on  Mach  number  and 
flow  angle  indicated  in  figure  4,  the  extent  of  the  range  of  circumferential  wave  numbers 
for  which  acoustic  response  disturbances  carry  energy  away  from  the  blade 
row  is  directly  proportional  to  the  excitation  frequency  and  inversely  proportional  to 
free-stream  speed. 

Curves  illustrating  the  behavior  of  the  axial  wave  number,  and  attenuation  con¬ 
stant,  /?,  for  various  subsonic  free-stream  Mach  numbers  and  a  free-stream  flow  angle  of 
45  deg  are  shown  in  figure  5.  Curves  for  M  =  0.8  and  flow  angles  of  15  deg,  45  deg  and 
75  deg  are  shown  in  figure  6.  The  points  on  these  curves  at  which  (5  fdt^  |M,n)  are  dis¬ 
continuous  are  the  resonant  values  of  £r?.  For  superresonant  disturbances, 
the  far  upstream  and  far  downstream  axial  wave  numbers  differ  and  0  =  0.  The  differ¬ 
ence  1^}00  increases  with  increasing  Mach  number.  For  subresonant  disturbances, 

°°?  £^)  or  the  far  upstream  and  far  downstream  wave  numbers  are 

identical  and  depend  linearly  on  £f?.  The  attenuation  constant  0  increases  with  \£n\ 
and  depends  linearly  on  \lft  |  as  l^j  — >  oo.  Thus,  except  at  Mach  numbers  close  to 
one,  acoustic  response  disturbances  of  high  wave-number  magnitude  attenuate  rapidly 
with  increasing  distance  from  the  blade  row.  This  is  an  important  feature  which  allows 
the  successful  application  of  numerical  field  methods  to  resolve  unsteady  flows  through 
cascades  operating  at  subsonic  inlet  and  exit  conditions. 

To  examine  a  typical  situation  in  more  detail  we  again  consider  the  case  uj  =  1, 
G  =  1,  V  —  1,  and  H  =  45  deg.  In  addition,  we  consider  an  unsteady  excitation  at 
o  —  0.5.  Such  an  excitation  produces  an  (n  —  0)  acoustic  response  disturbance  at 
M  —  0.7,  0.8  and  0.9  which  carries  energy  away  from  the  blade  row  in  both  the  far 
upstream  and  far  downstream  regions.  The  calculations  used  to  determine  the  curves 
in  figures  5  and  6  indicate  that  the  axial  wave  numbers  of  the  upstream  propagating 
disturbance  are  1.736,  2.373  and  3.245,  respectively,  and  those  for  the  downstream 
traveling  disturbance  are  —0.494,  —0.572  and  -0.639,  respectively.  The  n  —  ±1,  ±2, . . . 
acoustic  response  fluctuations  at  =  o  +  27rn  attenuate  with  distance  from  the  blade 
row  and,  except  for  the  n  =  1  disturbance  at  M  =  0.9,  they  attenuate  quite  rapidly. 
The  attenuation  constant  and  axial  wave  number  for  the  n  =  1  disturbance  (i.e., 
£f?  =  kn  i  =  o  +  2n  =  6.783)  at  M  =  0.9  are  0  —  0.662  and  k^TOO  =  =  5.580, 

respectively. 

In  numerical  calculations  of  the  unsteady  flow  through  a  cascade,  propagating  acous¬ 
tic  disturbances  and  those  which  attenuate  gradually  with  increasing  distance  from  the 
blade  row  must  be  accurately  represented.  Thus  for  our  example  configuration,  with 
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M  =  0.9,  acoustic  waves  with  k n  =  6.783  and  k ^  =  5.580  must  be  resolved  accu¬ 
rately.  In  this  case,  the  computational  mesh  requirements  are  not  stringent,  but  they 
can  become  so  if  either  the  Mach  number  or  the  excitation  frequency  is  increased.  In 
particular,  in  the  foregoing  example  if  M  =  0.95  and  a  —  Odeg,  propagating  (/ 3  =  0) 
acoustic  response  disturbances  with  knto  =  0,  knt\  =  6.28  and  kn)2  —  13.57  occur.  The 
axial  wave  numbers  of  these  propagating  disturbances  in  the  far  upstream  region,  k$  _oo, 
are  2.89,  10.43  and  13.76,  respectively.  Consequently,  small  mesh  spacings  are  required 
to  capture  the  n  —  1  and  n  =  2  acoustic  response  waves. 

5.2.2  Supersonic  Flow 

Resonance  curves  for  supersonic  free-stream  flows  with  subsonic  axial  velocity  com¬ 
ponent  are  shown  in  figure  7.  As  M  cos  fi  — >  1  from  below,  both  and  approach 
the  value  -M / (M2  —  l)1  2.  The  locus  of  such  points  is  shown  as  a  dashed  curve  in  figure 
7.  If  M  cos  VI  >  1,  the  axial  velocity  is  supersonic,  resonance  does  not  occur  and  acoustic 
response  disturbances  propagate  only  in  the  downstream  axial  direction.  For  super¬ 
sonic  flows  with  subsonic  axial  velocity  component,  the  resonant  condition  tn  = 
varies  substantially  with  Mach  number,  particularly  at  low  supersonic  Mach  numbers. 
Indeed,  — *  —  oo  as  M  — +  1  from  above  but,  as  in  subsonic  flow,  )  — »  —  (2sinn)  1 . 

The  extent  of  the  superresonant  region,  i.e.,  {^+\ decreases  with  increasing  Mach 
number  and  increases  with  increasing  flow  angle. 

As  a  specific  example,  consider  the  case  u>  =  1,  G  =  1,  V  =  1  and  0  =  60  deg.  The 
calculations  leading  to  the  curve  for  0  =  60  deg  indicate  that  for  a  free-stream  Mach 
number  of  1.6  and  interblade  phase  angles  in  the  range  [  —  7r,  7r]  all  acoustic  response 
disturbances  persist  in  the  far  field  and  carry  energy  away  from  the  blade  row  if  a  < 
—2.037  or  a  >  —0.806,  but  one  such  wave  will  attenuate  with  increasing  distance  from 
the  blade  row  (the  n  =  0  wave)  if  —2.037  <  o  <  —0.806.  For  a  Mach  number  of  1.4 
the  corresponding  ranges  are  a  <  —2.810,  a  >  —0.727  and  —2.810  <  a  <  —0.727. 
Finally,  for  a  Mach  number  of  1.2  the  n  =  0  acoustic  response  disturbance  will  be  of 
decaying  type  if  —  n  <  a  <  —0.652,  the  n  =  —  1  disturbance  will  be  of  decaying  type 
if  1.267  <  o  <  7r,  but  all  disturbances  are  of  propagating  type  if  —0.652  <  a  <  1.267. 
As  the  Mach  number  is  reduced  further,  additional  acoustic  disturbances  will  be  of 
decaying  type. 

Axial  wave  numbers  and  attenuation  constants  of  the  far-field  acoustic  response 
disturbances  for  Mach  numbers  of  1.1,  1.2,  1.3,  1.4  and  1.5  and  a  flow  angle  of  60deg 
are  shown  in  figure  8,  and  those  for  a  Mach  number  of  1.3  and  flow  angles  of  45  deg, 
60  deg  and  75  deg,  in  figure  9.  Acoustic  response  disturbances  at  subresonant  tangential 
wave  numbers,  i.e.,  £^e(— oo,  or  £f7e(£^  ^,oo),  persist  and  carry  energy  away  from 
the  blade  row.  Thus  acoustic  response  disturbances  of  high  wave-number  magnitude 
persist  in  the  far  field  and,  as  such,  they  can  impose  serious  limitations  on  the  numerical 
field  methods  developed  for  predicting  unsteady  flows  through  cascades  operating  at 
supersonic  inlet  and/or  exit  conditions. 
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5.3  The  Potential  Response  to  Vortical  Disturbances 

The  axial  and  tangential  wave  numbers  of  the  potential  response  in  the  far  down¬ 
stream  region  caused  by  continuous  rotational  velocity  fluctuations  can  also  have  high 
wave-number  content(see  (4.4)),  but  as  indicated  in  (4.22),  the  amplitude  of  such  po¬ 
tential  disturbances  varies  inversely  with  wave-number  magnitude.  The  potential  re¬ 
sponse  to  shed  (wake)  vorticity  has  a  wave-number  vector  of  —uV^er]  hence,  axial 
and  tangential  wave  numbers  of  k %  —  — wV^1  cos  fioo  and  kn  —  — sin  Hoq.  The 
magnitude  of  this  response  fluctuation  depends  upon  the  discontinuities  [<£]/?e/  and 
at  the  reference  wake  location  (£+,?7+).  The  potential  response  due  to  this  vor¬ 
ticity  has  a  rather  complicated  behavior  in  the  cascade  circumferential  direction  that 
depends  on  the  blade  spacing,  the  downstream  free-stream  velocity  and  the  temporal 
frequency  and  interblade  phase  angle  of  the  unsteady  excitation.  For  example,  the 
function  4>R {t+,ri) /M Ref  is  plotted  vs.  r\  -  r]+  in  figure  10  for  a  cascade  with  <2=1, 
operating  in  a  downstream  free  stream  with  Vqo  =  1  and  0^  =  45  deg  and  subjected 
to  non-vortical  excitations,  e.g.,  blade  motions  or  incident  acoustic  waves,  at  o  =  Odeg 
and  at  various  frequencies. 

At  the  temporal  frequencies  of  interest  for  practical  applications,  the  potential  dis¬ 
tribution  <^(£+,1?)  in  the  cascade  circumferential  direction  does  not  usually  impose 
stringent  requirements  on  the  spacing  of  mesh  lines  parallel  to  the  mean  flow  direction. 
However  for  certain  parametric  combinations,  it  can  become  necessary  to  pack  such 
lines  near  the  mean  wake  locations  to  resolve  the  potential  fluctuation  caused  by  wake 
vorticity. 
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6.  NUMERICAL  EXAMPLES:  SUBSONIC  FLOW 


Our  purposes  in  this  section  are  to  demonstrate  the  results  that  can  be  achieved 
by  matching  finite-difference  near-field  and  analytic  far-field  solutions  for  the  linearized 
unsteady  potential  and  to  illustrate  important  effects  associated  with  the  unsteady  flow 
behavior  in  the  far  field  on  the  aerodynamic  response  at  a  moving  blade  surface.  We  will 
restrict  our  consideration  to  two  simple  subsonic  cascade  configurations — a  flat-plate 
cascade  and  a  cascade  of  5%  thick  flat-bottomed  double-circular-arc  (DCA)  airfoils. 
In  each  case  the  cascade  stagger  angle  0  is  45  deg,  the  blade  spacing  G  is  unity  and 
the  inlet  free-stream  Mach  number  is  0.8.  Also, the  blades  undergo  pitching  oscillations 
about  midchord,  i.e.,  r  =  ax  i?r(see  (2.1)),  where  Rp  =  X  —  0.5ex.  Such  two  dimensional 
blade  motions  model  torsional  vibrations  of  actual  rotor  blades. 

Theoretical  results  will  be  presented  for  the  complex  amplitudes  of  the  first  har¬ 
monic  unsteady  aerodynamic  moments, 

m  =  ij^PgRp  ■  dr  ,  (6.1) 

and  pressure  difference  distributions, 

Ap(x)  =  pB  ( X ,  y-)-pB  (z,  y+),  (6.2) 

acting  on  the  reference  blade  of  each  cascade.  Here  ps  is  the  complex  amplitude  of 
the  first-harmonic  pressure  acting  at  the  moving  reference  blade  surface,  the  subscripts 
refer  to  the  upper  (+)  and  lower  (— )  surfaces  of  the  blade  and  the  integration  is  taken 
over  the  mean  blade  surface  (B).  Results  will  also  be  presented  for  the  velocity  potential 
distribution  along  the  stagnation  streamline  that  coincides  with  the  upper  surface  of 
the  reference  blade. 

We  consider  pitching  (torsional)  motions  with  a  =  1  at  four  different  frequencies, 
u>  —  0.5,  1.0,  1.5,  and  2.0,  and  at  interblade  phase  angles  lying  in  the  range  [ — 7r,  tt] . 
Since  a  is  taken  to  be  a  real  quantity,  the  real  and  imaginary  components  of  a  refer¬ 
ence  blade  complex  response  parameter  represent  components  that  are  in-  and  out-of- 
phase,  respectively,  with  the  reference  blade  displacement.  The  stability  of  a  torsional 
blade  motion  depends  upon  the  sign  of  the  out-of-phase  moment.  If  /m{oc}  =  0  and 
Im{m}  >  0,  the  airstream  supplies  energy  to  the  blade  motion  and  that  motion  is 
unstable  according  to  linearized  theory. 

The  steady  flows  through  the  flat  plate  and  DCA  cascades  are  assumed  to  satisfy  a 
zero-load  or  Kutta  condition,  i.e., 

V  ■  dr |B_  =  -V  ■  dr\D+  (6.3) 

at  the  leading  and  trailing  edges  of  each  blade.  Therefore  the  inlet  free-stream  Mach 
number  is  the  only  far-field  quantity  needed  to  completely  specify  the  mean  potential 
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flow  through  the  cascade.  In  particular,  the  steady  flow  through  the  example  flat-plate 
cascade  is  simply  a  uniform  stream  with  M  =  0.8  and  fl  =  45  deg.  For  the  DCA 
cascade  operating  at  M_ =  0.8  the  (calculated)  inlet  flow  angle,  exit  Mach  number 
and  exit  flow  angle  are  49.4  deg,  0.62  and  43  deg,  respectively.  This  flow  is  entirely 
subsonic  with  a  maximum  Mach  number  of  0.941  occurring  at  a  point  on  the  suction 
surface  which  is  36.5%  of  blade  chord  downstream  from  the  leading  edge.  The  blade 
surface  Mach  number  distribution  is  shown  in  figure  3  of  Verdon  &:  Caspar  (1984).  The 
wave-number  behavior  of  the  far-field  acoustic  response  for  the  two  example  cascades 
is  shown  in  figure  11. 

Although  linearized  unsteady  flow  solutions  for  the  example  flat-plate  cascade  can 
be  determined  using  semi-analytic  surface-integral  methods  (e.g.,  see  Smith  1971),  the 
results  presented  herein  for  both  the  flat-plate  and  DCA  cascades  have  been  determined 
using  the  finite-difference  numerical  model  described  in  Verdon  &:  Caspar  (1982,  1984). 
Those  determined  for  the  flat-plate  cascade  are  in  very  good  agreement  with  Smith’s 
solutions.  The  finite-difference  calculations  were  performed  on  a  mesh  composed  of  ax¬ 
ial  lines  (£  =  constant)  which  are  parallel  to  the  blade  row  and  tangential  curves  which 
are  percentile  averages  of  the  upper  and  lower  boundaries  of  the  extended  blade-passage 
solution  domain.  This  mesh  extended  one  axial  chord  upstream  and  downstream  from 
the  blade  row  (i.e.,  —1  <  £/  cos  0  <  2).  For  the  most  part  uniform  mesh  spacings  were 
used  with  A£  =  0.03  cos  0  and  Arj  =  3%  of  the  distance  between  the  upper  and  lower 
boundaries,  but  mesh  lines  were  concentrated  near  blade  edges  and  near  the  upper  and 
lower  boundaries  (i.e.,  the  blades  and  their  wakes)  of  the  solution  domain.  The  uniform 
spacings  given  above  were  selected  so  that  the  unsteady  flows  at  the  highest  frequency 
considered,  u>  =  2.0,  could  be  resolved  accurately  over  the  entire  interblade  phase  angle 
range  [  —  ir,  7rj. 

6.1  Flat-Plate  Cascade 

Unsteady  response  predictions  for  the  flat  plate  cascade  operating  in  a  uniform 
stream  with  M  =  0.8  and  fl  =  45  deg  are  shown  in  figures  12  through  15.  The 
resonant  circumferential  wave  numbers  in  both  the  far  upstream  and  far  downstream 
regions  for  this  configuration  are  k =  — 0.575cu  and  k W  —  3.09tu.  Acoustic  response 
disturbances  at  circumferential  wave  numbers,  knJl,  lying  between  these  values  persist 
in  both  far-field  regions  and  carry  energy  away  from  the  blade  row.  Disturbances  for 
which  kn  n  <  or  kn  n  >  kjf)  attenuate  with  increasing  distance  from  the  blade  row. 

Unsteady  moments  due  to  torsional  vibrations  about  midchord  at  c v  =  0.5,  1.0,  1.5 
and  2.0  are  shown  in  figure  12.  The  arrows  above  each  curve  indicate  the  resonant 
values  of  a  i.e.,  ffW  =  k^/G  —  2nn  where  n  =  0  for  \k^/G\  <  n,  n  =  1  for 
7r  <  k^/G  <  37T,  etc.,  and  therefore  the  boundaries  between  the  different  types  of 
unsteady  excitation.  Note  the  abrupt  changes  in  the  unsteady  moment  that  occur  near 
resonance  and,  more  importantly,  the  different  character  of  the  moment  response  in  the 
different  regions  of  unsteady  excitation.  For  u>  =  0.5  and  cu  =  1.0  the  blade  motions 
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are  superresonant  (1,1)  with  the  n  =  0  acoustic  response  disturbance  persisting  far 
upstream  and  far  downstream  of  the  blade  row  for  values  of  o  lying  between  the  two 
resonant  values.  The  blade  motions  at  u  =  0.5  and  u  =  1.0  are  subresonant  for 
cre[  — 7T,  ed  ))  and  ae(at+),7r],  At  u)  —  1.5  the  blade  motions  are  superresonant  (1,1) 
for  cre[ — 7r, at+))  and  ae(at_),7r],  and  subresonant  for  ae(at+),at~)).  In  the  first  case 
the  n  =  1  acoustic  response  disturbance  is  of  propagating  type  and  in  the  second  the 
n  =  0  disturbance  is  of  propagating  type.  The  flat-plate  blade  motions  at  cu  =  2  are 
superresonant  at  all  interblade  phase  angles.  These  motions  are  superresonant  (1,1) 
for  ae(— 7r,at-))  and  ae(at+),7r)  and  superresonant  (2,2)  for  ac(at-), at+)).  In  the  first 
case  the  n  —  1  acoustic  response  disturbance  is  of  propagating  type,  in  the  second  the 
n  =  0  disturbance  is  of  propagating  type  and  in  the  third  both  the  n  =  0  and  the  n  =  1 
acoustic  response  disturbances  are  of  propagating  type.  The  foregoing  information  is 
summarized  in  Table  1. 

Table  1.  Unsteady  far- field  behavior  for  the  flat-plate  cascade 


Excitation  frequency 
and 

resonant  phase  angles 

Interblade 
phase  angle 
range 

Type  of 
blade  motion 

Propagating  Waves 
upstream/downstream 

u>  =  0.5 

cr(_)  =  -0.288 

ore)  =  1.545 

(  tt,  crt  — )) 

(at-),  at+)) 
(at+),  7 r] 

subresonant 

superresonant(l,l) 

subresonant 

none/none 

n=0/n=0 

none/none 

u>  =  1.0 
cr(_)  =  -0.575 
ct(+)  =  3.090 

[— 7r,at-)) 
(crt-),  a*4') 
(a*+),  7r] 

subresonant 

superresonant 

subresonant 

none/none 

n=0/n=0 

none/none 

w  =  1.5 
at-)  =  -0.853 
a(+)  =  -1.647 

[— 7r,  at+)) 
(at+),at-)) 

(a(-),  7r] 

superreson  ant  (1,1) 
subresonant 
superresonant  (1,1) 

n=l/n=l 
none/none 
n=0/n— 0 

co  =  2.0 
at")  =  -1.150 
at+)  =  -0.103 

[-7r,at-)) 
(at-),  at+)) 
(at+),  7r  ] 

superresonant  (1,1) 
superresonant  (2, 2) 
superresonant  (1,1) 

n=l/n  =  l 

n=0,l/n=0,l 

n=0/n=0 

Unsteady  pressure  difference  distributions  for  the  example  flat-plate  cascade  with 
blades  vibrating  at  w  =  0.5,  1.0,  1.5  and  2.0  are  shown  in  figure  13  for  a  =  —  1  and 
in  figure  14  for  a  =  2.  The  blade  motions  at  a  —  -1  are  subresonant  at  the  three 
lower  frequencies  and  superresonant  (2,2)  at  ui  —  2.  In  the  latter  case  the  propagating 
acoustic  response  waves  have  wave  numbers  [k^kn)  of  (1.778,  -1)  and  (5.931,  5.283) 
in  the  far  upstream  region  and  (—0.577,-1)  and  (1.738,  5.283)  in  the  far  downstream 
region.  The  blade  motion  at  o  =  2  and  cu  =  0.5  is  subresonant.  The  motions  at 
o  =  2  and  u>  =  1.0,  1.5,  and  2.0  are  superresonant  (1,1)  with  wave  numbers  of  (3.085, 
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2),  (4.362,  2)  and  (5.474,  2),  respectively,  in  the  far  upstream  region  and  (0.129,  2), 
(-0.483,  2)  and  (-0.930,  2),  respectively,  in  the  far  downstream  region.  The  pressure 
difference  curves  for  a  =  -1  differ  substantially  for  the  four  frequencies  considered.  One 
reason  for  this  is  that  the  motions  at  o  =  —  1  occur  close  to  the  resonance  condition 
cr(-)  =  -0.575w,  particularly  those  at  the  two  higher  frequencies.  In  contrast,  the  blade 
motions  at  o  —  2  lie  well  within  the  subresonant  region  for  to  —  0.5  and  well  within 
the  superresonant  region  for  to  =  1.0,  1.5  and  2.0.  Thus,  although  the  magnitudes 
of  the  unsteady  pressure  differences  (and  the  unsteady  moments)  tend  to  be  much 
greater  for  the  motion  at  to  =  0.5  than  for  those  at  to  =  1.0,  1.5  and  2.0,  the  pressure 
difference  curves  for  a  =  2  show  a  consistent  trend  with  increasing  frequency.  Note 
that  at  both  phase  angles,  o  =  -1  and  a  =  2,  the  pressure  difference  curves  show 
higher  wave-number  content  with  increasing  temporal  frequency. 

Unsteady  potential  distributions  on  the  first  tangential  line  of  the  computational 
mesh  used  for  the  flat-plate  calculations  (i.e.,  on  y  =  0+)  are  shown  in  figure  15  for 
blade  motions  occurring  at  o  =  —  1  and  to  =  1,  o  =  2  and  to  —  0.5,  o  =  2  and  to  —  1  and 
at  o  =  -1  and  to  =  2.  The  results  for  -1  <  £/  cos  0  <  2  have  been  determined  from  the 
finite  difference  calculation,  those  for  £/  cos  0  <  -1  and  £/  cos  0  >  2,  from  the  analytic 
far-held  solutions.  The  two  solutions  were  matched  at  £/  cos  0  =  — 1  and  £/  cos  0  =  2. 
The  potential  in  the  far  upstream  region,  <f>  —  <f>p ,  is  due  solely  to  the  acoustic  response 
to  the  blade  motion,  while  that  in  the  far  downstream  region,  <f>  =  (f)P  +  (f>B ,  depends 
also  upon  the  vorticity  shed  from  the  blades.  The  potential  fluctuation  <f)p  has  a  wave 
number  in  the  streamwise  (in  this  case  the  x— )  direction  of  —  to.  The  blade  motions 
at  a  =  -1  and  to  -  1  and  at  a  =  2  and  to  =  0.5  are  subresonant,  so  <f>p  — ►  0  as 
|f|  oo.  In  both  cases  the  potential  fluctuation  in  the  far  downstream  region  is  due 
essentially  to  wake  vorticity.  However,  as  seen  from  the  potential  distribution  in  the 
far  upstream  region,  the  n  —  0  acoustic  response  wave  for  the  blade  motion  at  o  2 
and  to  =  0.5  attenuates  only  gradually  with  increasing  distance  from  the  blade  row. 
The  acoustic  response  to  the  superresonant  (1,1)  blade  motion  at  o  —  2  and  to  —  1  has 
a  streamwise  wave  number  of  kx  =  3.596  in  the  far  upstream  region  and  kx  =  1.505 
in  the  far  downstream  region.  The  potential  fluctuation  due  to  the  acoustic  response 
in  the  far  downstream  region  occurs  at  a  much  smaller  amplitude  than  the  fluctuation 
due  to  wake  vorticity.  Finally,  the  acoustic  response  to  the  superresonant  (2,2)  blade 
motion  at  a  =  -1  and  to  =  2  has  streamwise  wave  numbers  of  0.550  and  7.92  in  the 
far  upstream  region  and  -1.115  and  4.965  in  the  far  downstream  region.  The  acoustic 
disturbances  at  the  higher  wave  numbers  occur  at  much  smaller  amplitude  than  those 
at  the  lower  wave  numbers. 

6.2  The  DCA  Cascade 

Unsteady  response  predictions  for  the  example  DCA  cascade  are  presented  in  fig¬ 
ures  16  through  19.  Since  the  inlet  (M_ «,  =  0.8,n_oo  =  49.4  deg)  and  exit  (MTO  = 
0.62,  Uoo  =  43  deg)  free-stream  conditions  differ  for  this  configuration,  the  acoustic 
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response  characteristics  in  the  far  upstream  and  far  downstream  regions  also  differ. 
In  particular,  the  resonant  circumferential  wave  numbers  are  k^l^  =  —  0.547a;  and 
=  3.247a;  in  the  far  upstream  region  and  =  -0.595a;  and  k^  =  1.669a; 
in  the  far  downstream  region.  Acoustic  response  disturbances  occurring  at  circumfer¬ 
ential  wave  numbers,  kniU,  lying  between  k^  and  0  persist  in  the  far  upstream 
region,  and  those  occurring  at  circumferential  wave  numbers  lying  between  and 

^rj+oo  persist  in  the  far  downstream  region. 

The  unsteady  moments  over  the  interblade  phase  angle  range  [— 7r,7r]  are  shown  in 
figure  16  for  torsional  blade  vibrations  at  w  =  0.5,  1.0,  1.5  and  2.0.  Again  the  arrows 
(now  four)  above  each  curve  indicate  the  resonant  interblade  phase  angles.  Because 
of  the  different  acoustic  environments  in  the  two  far-field  regions,  it  is  a  more  difficult 
task  to  classify  the  DCA  blade  motions.  For  example,  consider  the  blade  motions 
at  w  =  1.  Those  occurring  in  the  range  are  subresonant,  and  therefore 

the  acoustic  response  to  the  blade  motion  attenuates  with  increasing  distance  from 
the  blade  row.  Blade  motions  at  the  other  non-resonant  interblade  phase  angles  are 
superresonant.  In  particular,  these  motions  are  superresonant  (1,0),  (0,1),  (1,1)  and 
(1,0)  for  interblade  phase  angles  lying  in  the  ranges  [-7T, erL"t2),  (^oo^a-oo)»  (^-oL *) 
and  (<tW,7t],  respectively.  In  the  first  case  the  n  —  1  acoustic  response  disturbance 
persists  in  the  far  upstream  region;  in  the  second,  the  n  =  0  disturbance  persists  far 
downstream;  in  the  third,  the  n  =  0  disturbances  persist  in  both  far-field  regions;  and, 
finally,  in  the  fourth,  the  n  =  0  disturbance  persists  far  upstream.  This  information  is 
summarized  in  Table  2  below  along  with  the  corresponding  results  for  the  DCA  blade 
motions  at  w  =  0.5,  1.5  and  2.0.  There  are  significant  differences  between  the  unsteady 
moment  responses  to  the  fiat-plate  and  DCA  blade  motions.  These  differences  can  be 
attributed  to  the  effects  of  mean  flow  variations  near  the  surfaces  of  the  DCA  blades 
and  to  the  different  far-field  acoustic  response  environments  seen  by  the  two  cascades. 

Unsteady  pressure  difference  distributions  for  the  DCA  cascade  are  shown  in  figures 
17  and  18  for  torsional  blade  motions  at  o  =  —1  and  o  —  2,  respectively.  The  motions 
at  a  =  -1  are  subresonant  for  u  =  0.5,  1.0  and  1.5  and  superresonant  (2,1)  for  u>  =  2.0. 
In  the  superresonant  case  the  wave  numbers  {k^kn)  of  the  propagating  acoustic  distur¬ 
bances  are  (1.403,-1)  and  (5.720,  5.283)  in  the  far  upstream  region  and  (-0.245,-1) 
in  the  far  downstream  region.  For  <7  =  2  the  blade  motion  at  a;  =  0.5  is  subresonant 
and  those  at  oj  =  1.0,  1.5  and  2.0  are  superresonant.  The  superresonant  (1,0)  motion 
at  <j  =  2  and  to  =  1.0  produces  an  upstream  propagating  acoustic  response  disturbance 
at  (kt,kn)  =  (2.906,2).  For  the  superresonant  (1,1)  motion  at  <7  =  2  and  u>  =  1.5  the 
wave  numbers  of  the  acoustic  response  are  (4.067,  2)  far  upstream  and  (-0.039,2)  far 
downstream.  Finally,  for  the  superresonant  (1,1)  blade  motion  at  cr  =  2  and  u>  =  2  the 
propagating  acoustic  disturbances  occur  at  (5.080,  2)  far  upstream  and  (—0.665,2)  far 
downstream. 
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Table  2.  Unsteady  far-field  behavior  for  the  DCA  cascade 
Excitation  frequency  Interblade  Type  of  Propagating  Waves 

and  phase  angle  blade  motion  upstream/downstream 

resonant  phase  angles  range  


(jj  = 

0.5 

[“*><22) 

subresonant 

none/none 

W  —  OO 

=  -0.274 

(<22><22) 

superresonant  (0,1) 

none/n— 0 

=  1.624 

(<22,<22) 

superresonant  (1,1) 

n=0/n=0 

=  -0.297 

(<22><22) 

superresonant  ( 1 ,0) 

n— 0/n=0 

a(+) 

OO 

=  0.834 

(<22,2 

subresonant 

none/none 

0J  = 

1 

(-7T,<22) 

superresonant(l,0) 

n~]  /none 

2 

=  -0.547 

(<22,<22) 

subresonant 

none/none 

<22 

=  -3.037 

(<22,<22) 

superresonant  (0,1) 

none/n— 0 

(j(~) 

OO 

=  -0.595 

(22,<22) 

superresonant(l,l ) 

n— 0/n=0 

<j(  +  ) 

OO 

=  1.669 

(<22>2 

superresonant  ( 1 ,0) 

n— 0/none 

U)  = 

1.5 

[-^,<22) 

superresonant  (1,0) 

n=l/none 

<22 

=  -0.821 

(<22,<22) 

subresonant 

none/none 

<22 

=  -1.413 

(<22><22) 

superresonant  (0,1) 

none/n— 0 

<22 

=  -0.892 

(<22,<22) 

superresonant  (1,1) 

n— 0/n— 0 

a(+) 

OO 

=  2.503 

(<22,2 

superresonant  (1,0) 

n— 0/none 

UJ  = 

2 

[-*■,<22) 

superresonant  (1,1) 

n=l/n=l 

<22 

=  -1.095 

(<22,<22) 

superresonant(l,0) 

n=l/none 

<22 

=  0.218 

(<22,<22) 

superresonant  (1,1) 

n  =  l/n=0 

OO 

=  -1.190 

(<22,<22) 

superresonant  (2,1) 

n— 0,1/n— 0 

cr(+) 

OO 

=  -2.947 

(<22,2 

superresonant  (1,1) 

3 

II 

o 

2 

II 

o 

The  pressure  difference  curves  for  the  DCA  and  flat-plate  blades  vibrating  at  a  =  —  1 
show  similar  behaviors  for  u>  =  0.5  and,  except  near  the  blade  leading  edge,  w  =  1.0. 
Here,  the  differences  between  the  DCA  and  flat-plate  results  are  due  primarily  to  the 
nonuniform  mean  flow  over  the  DCA  blade  surfaces.  The  responses  to  the  DCA  and 
flat-plate  blade  motions  at  a  =  —  1  differ  considerably  for  w  =  1.5  and  w  =  2.  These 
motions  occur  near  a  resonance  condition  and,  as  we  shall  see  more  clearly  below,  the 
flat-plate  and  DCA  blade  motions  produce  very  different  acoustic  responses  in  the  far 
field.  The  pressure  difference  curves  for  the  DCA  and  flat-plate  blade  motions  at  o  =  2 
are  quite  similar  for  the  subresonant  blade  motions  at  u;  =  0.5  but  differ  considerably 
for  the  superresonant  motions  at  w  =  1.0,  1.5  and  2.0.  Again,  the  differences  can  be 
attributed  to  the  nonuniform  mean  flow  over  the  surface  of  each  DCA  blade  and  to  the 
different  far  downstream  acoustic  response  environments  seen  by  the  two  cascades.  For 
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example,  at  a  =  2  and  u  =  1  an  acoustic  response  fluctuation  persists  far  downstream 
of  the  flat-plate  cascade,  but  all  such  disturbances  attenuate  with  increasing  distance 
downstream  of  the  DCA  cascade. 

Unsteady  potential  distributions  on  the  lower  boundary  of  the  extended  blade  pas¬ 
sage  solution  domain  (i.e.,  the  stagnation  streamline  that  coincides  with  the  upper 
surface  of  the  reference  blade)  are  shown  in  figure  19  for  DCA  blade  motions  at  o  =  —  1 
and  lj  —  1,  o  =  2  and  u  =  0.5,  a  —  2  and  to  =  1  and  for  o  =  — 1  and  w  =  1.  Again,  the 
results  for  —  1  <  £/  cos  0  <  2  have  been  determined  from  the  finite-difference  near-field 
solution  and  those  for  £/  cos  0  <  —  1  and  £/  cos  0  >  2  from  the  analytic  far-field  solu¬ 
tions.  The  far  downstream  potential  for  the  first  three  of  these  motions  is  essentially 
due  to  wake  vorticity  and  therefore  has  a  wave  number  of  =  —1.26a;  in  the 

streamwise  direction.  The  DCA  blade  motions  at  o  —  -1  and  w  =  1  and  at  a  —  2  and 
oj  =  0.5  are  subresonant,  and  the  far-field  behavior  associated  with  these  motions  is 
similar  to  that  observed  for  the  flat  plate  cascade.  The  far  upstream  acoustic  response 
to  the  superresonant  (1,0)  blade  motion  at  a  =  2  and  u>  —  1  occurs  at  a  streamwise  wave 
number  of  3.409  and  at  a  much  larger  amplitude  than  the  upstream  propagating  wave 
generated  by  the  corresponding  flat-plate  blade  motion.  Finally,  the  acoustic  response 
to  the  superresonant  (2,1)  motion  at  a  =  -1  and  u  =  2  has  streamwise  wave  num¬ 
bers  of  0.154  and  7.734  in  the  far  upstream  region  and  —0.861  in  the  far  downstream 
region.  As  can  be  seen  from  the  results  in  figures  15  and  19,  the  far-field  response 
associated  with  this  motion  differs  considerably  from  that  produced  by  the  fiat-plate 
blade  vibrations  at  o  =  —  1  and  u)  =  2. 
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7.  CONCLUDING  REMARKS 


Analytic  solutions  for  the  unsteady  velocity  potential  response  far  upstream  and 
far  downstream  of  an  isolated  two-dimensional  cascade  and  for  the  rotational  velocity 
response  far  downstream  have  been  provided.  Fourier  series  representations  account 
for  the  continuous  response  fluctuations  associated  with  acoustic  disturbances,  which 
either  attenuate  with  increasing  axial  distance  from  the  blade  row  or  persist  and  carry 
energy  away  from  or  parallel  to  the  blade  row,  and  the  rotational  velocity  distur¬ 
bances  which  are  convected  through  the  blade  row.  Closed  form  expressions  account 
for  the  discontinuous  response  fluctuations  associated  with  the  concentrated  vorticity 
convected  along  blade  wakes.  The  complete  far-field  solution  provides  the  upstream 
and  downstream  information  needed  for  formulating  the  boundary-value  problem  that 
describes  the  linearized  unsteady  perturbation  of  an  irrotational  and  isentropic  mean 
flow  through  a  cascade.  In  addition,  this  complete  solution  can  be  easily  implemented 
into  the  numerical  analyses  and  computer  codes  developed  to  predict  the  aerodynamic 
response  of  the  blading  to  prescribed  small-amplitude  unsteady  excitations,  i.e.,  blade 
motions  and  incident  entropic,  vortical  and  acoustic  disturbances. 

Special  cases  of  the  far-field  solutions  presented  herein  have  been  used  for  some 
time  in  unsteady  aerodynamic  analyses  intended  for  blade  flutter  prediction  in  which 
prescribed  blade  motions  are  the  only  source  of  unsteady  excitation.  However,  the 
solutions  provided  in  this  report  apply  more  generally  to  forced  vibration  problems  in 
which  blade  motions  caused  by  incident  external  aerodynamic  disturbances  are  also  of 
concern.  The  Fourier  series  representation  of  the  far-field  acoustic  response  is  particu¬ 
larly  well-suited  for  cascades  operating  at  subsonic  inlet  and  exit  Mach  number  since 
only  a  limited  number  of  acoustic  response  disturbances  can  persist  in  the  far  field. 
However,  as  indicated  by  our  numerical  results  for  flat-plate  and  DCA  cascades,  there 
is  a  complex  and  interesting  variety  of  far-field  response  phenomena  even  for  subsonic 
flows.  This  phenomena  must  be  represented  accurately  in  order  to  predict  correctly 
the  aerodynamic  response  at  a  vibrating  blade  surface.  At  high  subsonic  inlet  and  exit 
Mach  numbers  and  at  the  high  excitation  frequencies  typical  of  forced  blade  vibrations, 
stringent  mesh  spacing  requirements  must  be  imposed  in  numerical  calculations  of  un¬ 
steady  flows  through  blade  rows.  For  a  supersonic  free-stream  flow  with  subsonic  axial 
velocity  component  an  infinite  number  of  acoustic  response  disturbances  will  persist 
in  the  far  field.  Since  only  a  finite  number  of  these  disturbances  can  be  modeled  on  a 
computational  grid,  the  usefulness  of  the  Fourier  series  representation  of  the  far-field 
acoustic  response  for  such  flows  must  still  be  carefully  evaluated. 
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Figure  1.  Two-dimensional  subsonic  compressor  cascade;  <  M_oo  <1 
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Figure  2.  Far-field  acoustic  response  behavior  for  subsonic  flow;  the  first  number 
listed  under  each  region  refers  to  the  inequality  satisfied  by  k^p,  the  second  to 

the  expression  used  to  determine  Xn- 
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Figure  3.  Far-field  acoustic  response  behavior  for  supersonic  flow;  the  first  number 
listed  under  each  region  refers  to  the  inequality  satisified  by  k_n,  the  second  to 

the  expression  used  to  determine  xn- 
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Figure  4.  Resonance  conditions  for  an  isolated  blade  row  operating  in 

a  subsonic  free  stream. 
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Figure  5.  Effect  of  Mach  number  on  the  far-field  acoustic  response  of  an 
isolated  blade  row  operating  in  a  subsonic  free  stream;  Q  =  45  deg: 

(a)  axial  wave  number  vs.  circumferential  wave  number;  (b)  attenuation 
constant  vs.  circumferential  wave  number. 
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Figure  6.  Effect  of  flow  angle  on  the  far-field  acoustic  response  of  an 
isolated  blade  row  operating  in  a  subsonic  free  stream;  M  =  0.8:  (a)  and 

(b)  as  in  figure  5. 
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Figure  7.  Resonance  conditions  for  an  isolated  blade  row  operating  in 

a  supersonic  free  stream. 
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.  Effect  of  Mach  number  on  the  far-field  acoustic  response  of 
ted  blade  row  operating  in  a  supersonic  free  stream;  Q  =  60  deg: 
wave  number  vs.  circumferential  wave  number;  (b)  attenuation 
constant  vs.  circumferential  wave  number. 
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Figure  9.  Effect  of  flow  angle  on  the  f ar-f ielcf  acoustic  response  of 
an  isolated  blade  row  operating  in  a  supersonic  free  stream; 

M  =  1 .3:  (a)  and  (b)  as  in  figure  8. 
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Figure  10.  The  potential  response  to  concentrated  wake  vorticity; 
G  =  1,  Vro  =  1,  Qoo  =  45  deg,  o  =  0. 
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Figure  11.  Axial  vs.  circumferential  wave  number  of  the  far-field 
acoustic  response  associated  with  the  example  cascades;  G  =  1 , 
0  =  45  deg:  (a)  flat-plate  cascade;  M  =  0.8,  Q  =  45  deg;  (b)  DCA 
cascade;  M  _  oo  =  0.8,  Q_oo  =  49.4  deg,  Moo  =  0.62,  QTC  =  43  deg. 
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Figure  12.  Unsteady  moment  vs  interblade  phase  angle  for  torsional 

blade  vibrations  of  the  example  flat-plate  cascade: - in-phase 

component  (real  part);  - out-of-phase  component  (imaginary  part). 
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Figure  12.  (cont.)  Unsteady  moment  vs  interblade  phase  angle  for 

torsional  blade  vibrations  of  the  example  flat-plate  cascade: - 

in-phase  component  (real  part); - out-of-phase  component 

(imaginary  part). 
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Figure  13.  Unsteady  pressure  difference  distributions  due  to  torsional 
blade  vibrations  at  o  =  - 1  for  the  example  flat-plate  cascade:  (a)  in-phase 
component;  (b)  out-of-phase  component. 
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Figure  14.  Unsteady  pressure  difference  distributions  due  to  torsional 
blade  vibrations  at  o  =  2  for  the  example  flat-plate  cascade:  (a)  in-phase 
component;  (b)  out-of-phase  component. 
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Figure  15.  Unsteady  potential  distributions  on  y  =  0+  for  torsional 
blade  vibrations  of  the  example  flat-plate  cascade;  numerical 

solution  domain  extends  over  -1  <  |/cos  0  <2: - in-phase 

component;  - out-of-phase  component. 


87-3-27-17 


R87-957333-1 


Figure  15.  (cont.)  Unsteady  potential  distributions  on  y  =  0+  for 
torsional  blade  vibrations  of  the  example  flat-plate  cascade; 
numerical  solution  domain  extends  over  -1  <|/cos  ©  <2: 
- in-phase  component;  - out-of-phase  component. 
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Figure  16.  Unsteady  moment  vs  interblade  phase  angle  for  torsional 

blade  vibrations  of  the  example  DCA  cascade: - in-phase 

component  (real  part); - out-of-phase  component  (imaginary  part). 
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Figure  16.  (cont.)  Unsteady  moment  vs  interblade  phase  angle  for 

torsional  blade  vibrations  of  the  example  DCA  cascade: - 

in-phase  component  (real  part);  - out-of-phase  component 

(imaginary  part). 
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Figure  17.  Unsteady  pressure  difference  distributions  due  to  torsional 
blade  vibrations  at  o=  -1  for  the  example  DCA  cascade:  (a)  in-phase 
component;  (b)  out-of-phase  component. 


87-3-27-21 


Im  {Ap}  Re  I  AP  S 


R87-957333-1 


5 


(a) 


Figure  18.  Unsteady  pressure  difference  distributions  due  to  torsional 
blade  vibrations  at  o  =  2  for  the  example  DCA  cascade:  (a)  in-phase 
component;  (b)  out-of-phase  component. 
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Figure  19.  Unsteady  potential  distributions  along  the  reference 
stagnation  streamline  for  torsional  blade  vibrations  of  the 
example  DCA  cascade;  numerical  solution  domain  extends 

over  -1<  4/cos  0<  2: - in-phase  component; - 

out-of-phase  component. 
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Figure  19.  (cont.)  Unsteady  potential  distributions  along  the 
reference  stagnation  streamline  for  torsional  blade 
vibrations  of  the  example  DCA  cascade;  numerical  solution 

domain  extends  over  - 1  <  |/cos  0  <  2: - in-phase 

compoent;  _ out-of-phase  component. 
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